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Abstract. In this paper we investigate the evolution of a Jordan-Brans-Dicke scalar field, <h, 
with a power-law potential in the presence of a second scalar held, (f>, with an exponential 
potential, in both the Jordan and the Einstein frames. We present the relation of our model 
with the induced gravity model with power-law potential and the integrability of this kind 
of models is discussed when the quintessence held <j) is massless, and has a small velocity. 
The fact that for some fine-tuned values of the parameters we may get some integrable 
cosmological models, makes our choice of potentials very interesting. We prove that in 
Jordan-Brans-Dicke theory, the de Sitter solution is not a natural attractor. Instead, we 
show that the attractor in the Jordan frame corresponds to an “intermediate accelerated” 
solution of the form a(t ) ~ e aitn , as t —>• oo where a.\ > 0 and 0 < pi < 1, for a wide 
range of parameters. Furthermore, when we work in the Einstein frame we get that the 
attractor is also an “intermediate accelerated” solution of the form a(t) ~ e“ 2tP2 as t —>• oo 
where «2 > 0 and 0 < p 2 < 1, for the same conditions on the parameter space as in the 
Jordan frame. In the special case of a quadratic potential in the Jordan frame, or for a 
constant potential in the Einstein’s frame, the above intermediate solutions are of saddle 
type. These results were proved using the center manifold theorem, which is not based 
on linear approximation. Finally, we present a specific elaboration of our extension of the 
induced gravity model in the Jordan frame, which corresponds to a particular choice of a 
linear potential of <h. The dynamical system is then reduced to a two dimensional one, and 
the late-time attractor is linked with the exact solution found for the induced gravity model. 
In this example the “intermediate accelerated” solution does not exist, and the attractor 
solution has an asymptotic de Sitter-like evolution law for the scale factor. Apart from some 
fine-tuned examples such as the linear, and quadratic potential U($>) in the Jordan frame, it 
is true that “intermediate accelerated” solutions are generic late-time attractors in a modified 
Jordan-Brans-Dicke theory. 
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1 Introduction 

A large amount of research has been devoted to the explanation of the late-time accelera¬ 
tion of the universe, either by introducing the concept of Dark Energy or by modifying the 
gravitational sector itself. Among the simplest candidates for Dark Energy, one can find 
canonical scalar fields, phantom Helds or the combination of both Helds in a unified quintom 
model, see [1, 2]; for the second approach there are several attempts reviewed in [3] (see 
references therein). Despite their interpretation, both approaches can be transformed one 
into the other, since the crucial issue is just the number of degrees of freedom beyond General 
Relativity and standard model particles (see [4] for a review on such a unified point of view). 
Finally, the above scenarios are well-suited not just for late-time implications, likewise for 
the description of an inflationary stage [5]. 

One example of modified gravitational theory is the so called scalar-tensor theory of 
gravity [6-8], in particular the Jordan-Brans-Dicke (JBD) theory [6, 7]. In this theory the 
effective gravitational coupling is time-dependent. The strength of this coupling is determined 
by a scalar field, the so-called JBD field, oc G . In modern context, JBD theory appears 
naturally in supergravity models, Kaluza-Klein theories and in all known effective string 
actions [9-16]. Furthermore, we can promote the Brans-Dicke (BD) parameter, u>o, presents 
in the original theory to a non-constant BD parameter u;o(4>), and to consider a non-zero 
self-interaction potential £/(<&), even surviving astrophysical tests [17, 18]. 

In [19] it was investigated the dynamics of the JBD scalar field with a quadratic poten¬ 
tial and barotropic matter. The authors used the dynamical systems approach, revealing that 
the complexity of dynamical evolution, in homogeneous and isotropic cosmological models, 
depends on the BD parameter, u>q, and the barotropic matter index, w m . The authors claim 
that the quadratic potential function leads naturally to a de Sitter state. The results in [19] 
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were extended by [20] for an arbitrary potential function. In [21], it was investigated the ob¬ 
servational constraints on the JBD cosmological model using observational data coming from 
distant supernovae type la, the Hubble function, H(z) measurements, information coming 
from the Alcock-Paczyriski test, and baryon acoustic oscillations. However, the values found 
in [19-21] for the BD parameter ujq are several orders of magnitude lower than the bound 
cao > 4 x 10 4 imposed by the Solar System tests [22, 23], and the bounds estimated on the 
basis of cosmological arguments wo > 120 [24] and 10 < wo < 10 7 [25]. This was the main 
objection to the models [19-21] in [26]. 

In this latter paper, the authors stated that the de Sitter solution is an attractor in the 
Jordan frame of the BD theory only for the quadratic potential f7(<h) oc <b 2 . This result lead 
them to the claim that de BD cosmology does not have the ACDM model as the universal 
attractor. Additionally, the authors showed that in the stable de Sitter critical point, as 
well as in the stiff-matter equilibrium configurations, the dilaton is necessarily massless. Due 
to the recent discussions in the literature concerning this topic, we consider it is worthy 
to investigate the subject further. In this paper we investigate a JBD scalar field, 4>, with 

potential t/(<J>) = Uq& , 7 _i = i/wo + | (where we have chosen the positive square 
root by convention), in the presence of a second scalar held, <j>, with exponential potential 
V(<f)) = Voe~ Xv ‘t > to the matter content. Because the addition of 0, we call this scenario a 
“modified JBD” theory. We assume that the BD parameter wo is finite, so the limiting case 
wo —>• Too and <b —» 1 (we use 8itG = 1), where GR is recovered, is excluded here. 

Several gravity theories consider multiple scalar fields, e.g., assisted inflation scenarios 
[27-32], quintom dark energy paradigm [33-43], among others [44-55]. In particular there 
are some theories where the role of dark matter is played by a scalar held which dynamically 
behaves as dust during certain epoch in evolution [56-58]. The main motivation of this 
work is to analyse if the de Sitter solution represents a natural attractor in the modified 
JBD model presented. We investigate the Jordan and Einstein frames, and in both cases 
we prove that, under the parameter choices wo > — | and A jj < 0, the late time attractor 
has an effective equation of state parameter u>tot = — 1. This region in the parameter space 
is compatible with the ranges described by observations [22-25, 59]. We prove that in this 
modified JBD model, the de Sitter solution is not a natural attractor. Instead, we show that 
the attractor in the Jordan frame corresponds to an “intermediate accelerated” solution of 
the form a(t) ~ e° ltP1 as t —>• oo with a\ > 0, 0 < p\ < 1. Furthermore, when we work in the 
Einstein frame we get that the attractor is, as well, an “intermediate accelerated” solution 
of the form a(t) ~ e“ 2tP2 as t —>• oo with «2 > 0, 0 < p 2 < 1. 

A scale factor of the form a(t) = exp (At-f) where A > 0 and 0 < / < 1 was introduced 
in [60-62] in the context of inflation. Since the expansion of the universe with this scale 
factor is slower than the de Sitter inflation (a(t) = exp (Ht) where H is constant), but 
faster than the power-law inflation (a(t) = t q where q > 1), it was called intermediate 
inflation. Intermediate inflationary models arise in the standard inflationary framework as 
exact cosmological solutions in the slow-roll approximation to potentials that decay with 
inverse power-law of the inflaton field [63]. These models have been studied in some warm 
inflationary scenarios [64-73]. 

Intermediate inflation is also found in the context of scalar tensor theories with a variable 
BD field in the Jordan frame and different matter content. In Ref. [74] a fluid with constant 
state parameter is considered. In Ref. [75] the author takes into account a scalar field as 
matter source and intermediate inflation is found in the slow-roll approximation. In the 
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reference [67, 68] it was investigated warm intermediate inflation in the JBD theory but 
formulated in the Einstein frame. Since this kind of solutions appear as late-time attractors 
in our context, we call them “intermediate accelerated” solutions. 

We note that under the scalar field rescaling a = 2y/u)o& for wo > 0 and without 
the second scalar field cj), we obtain from our model a special case of the so-called induced 
gravity model, which is integrable for a power-law potential £/(<£). The general solution of 
this class of models is known, see for example [76]. After conformal transformation, the 
induced gravity model with power-law potential becomes a General Relativity model with 
an exponential potential, that it is integrable as well [76, 77]. In this sense, our model can 
be considered a generalization of the induced gravity models described above since we have 
included a new scalar field < f> as the matter source. So it would be interesting to see how 
the behavior of the solutions for the induced gravity model changes when a small scalar field 
is added. In the subsection 2.1 a discussion of this issue is presented. However, the main 
purpose of our investigation is not to find analytical solutions but to study the asymptotic 
behavior of the solutions space of this kind of scenarios without using fine-tuning of the 
parameters and initial conditions. Dynamical systems theory is a powerful tool for doing this 
research. Nevertheless, the fact that for some fine-tuned values of the parameters we may 
get some integrable cosmological models, makes our choice of potentials very interesting. 

The paper is organized as follows. In section 2 the model is presented and the field 
equations in the Jordan’s frame are displayed. In the subsection 2.1 the relation of our model 
with the induced gravity model for power-law potential is presented and the integrability of 
this kind of models is discussed when the quintessence field cj) is massless and has a small 
velocity, <f>. In the subsection 2.2, the system is written as a dynamical system and the 
stability of the critical points is discussed. We separated the analysis in two parts: the analysis 
at the finite region, and the analysis at the infinite region, covering all the possibilities. 
In the subsection 2.3 we present the intermediate accelerated solution as a possible future 
attractor in the Jordan frame. In subsection 2.4 we investigate a linear potential of <h, 
which is equivalent to an extension of the so-called induced gravity model [76, 78]. The 
dynamical system is reduced to a two dimensional one, and the late-time attractor is linked 
with the solutions found in section 2.1. In section 3 the equations are written in the Einstein’s 
frame, through a conformal transformation. In the subsection 3.1 the system is written as 
a dynamical system and the stability of the critical points is discussed. Special emphasis 
is given to the possible late time attractors. In subsection 3.2 we present the intermediate 
accelerated solution as a possible future attractor in the Einstein frame. Concluding remarks 
are given in section 4. 

2 Field Equations in the Jordan’s frame 

Let us consider the action written in the Jordan frame as given by: 

SjF = f V=g - U (4>) - l -g^d^d„ct> - V(cj>)\ d 4 x, (2.1) 

where $ denotes the JBD scalar field, ujq is the BD parameter and cj) represents a quintessence 

2 —hlL 

scalar field. For the sake of simplicity we restrict our attention to the cases U($>) = Uo& ~i , 
with 7 -1 = \Jojq T ^ and V(4>) = Voe _Av< ^, but the analysis can be extended to general 
potentials using similar techniques as in [19, 26]. A jj and Ay are constants. By construction 
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we have assumed 7 is positive and finite, it follows uq > — | (the value uq = — | gives 7 
infinity). The JBD scalar field plays the role of an effective Planck mass, consequently we 
assume > 0. However, it can asymptotically evolves to its minimum value T = 0 as we 
shall show in the following sections. 

By considering a flat Friedmann-Lemaitre-Robertson-Walker (FLRW) metric: 


ds 2 = — dt 2 + a(t) 2 [dr 2 + r 2 (d,0 2 + sin 2 9dip 2 )] , 
the field equations become 

_ ( 12 ^ 12H 2 d> 24>H'(<1>) 2a;o<h 2 3 </> 2 

y2(Uo + 3 ) 2wo + 3 2loq + 3 (2wo + 3)4? 2ujq + 3 

4> = -3H<j>-V'((/>), 

3 H 2 <f> = + [/($) + 1^2 + _ 3 #^ 

^ 6 H 2 t/ 7 ^) wo( 2 u;o + 1)4 > 2 wo ^> 2 

(2a ; 0 + 3)<f> 2a ; 0 + 3 + 2w 0 + 3 2(2u ; 0 + 3)<F 2 (2u ; 0 + 3)<F’ 


( 2 . 2 ) 


(2.3a) 

(2.3b) 

(2.3c) 

(2.3d) 


where the Hubble expansion rate is given by H = ^ and the dot denotes derivatives with 
respect to the cosmic time. 

By defining the following effective energy densities p \, P 2 and the effective pressuress 
Pl,P2- 


pi = 3H 2 (1 - 4>) + + [/($) - 3 H$, 

P2 = \<t > 2 + V(<j)), 

( 8w 0 \ ^ , H 2 ((6w 0 + 9) <F - 6w 0 + 3) 

P\ = H \ 6 — --—— 4> H- 

\ (2wo + 3) + y 2cuo + 3 


/ 2cUq + cjq _ wo_\ ^2 1 f 2o;q 
V (2cj 0 + 3) <F 2 2$/ \ (2cj 0 + 3) <F 

P 2 = \(t> 2 -v{4>), 


the system (2.3) 


can be written as 



(2.4a) 


2 t/'(<F) 

2loq + 3 


(2.4b) 


U($)+ 


(2.4c) 


(2.4d) 


pi + 3H(pi +pi) = 0, (2.5a) 

P 2 + 3H(p 2 +P 2 ) = 0, (2.5b) 

n 2 = \ (/+ + P*), (2.5c) 

H =-^(pi+pi+P 2 +P 2 ) ■ (2.5d) 


The above phenomenological definitions of the energy densities are not unique, specially if 
an interaction term between both fields is considered [79]. 
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2.1 Relation with the induced gravity model 

Let us observe that by setting p$ := ^<j) + V{4>) = 0 , we obtain from ( 2 . 1 ) the so-called 
induced gravity model [76, 78]: 

SiG = J (w (cr)R - d^advCT - U (cr 2 )^ d A x, (2.6) 

under the choices a = 2y/ujo& and W(a) = given uj q > 0. This model admits exact 

solutions that we want to discuss in the following. 

2—ha i / r 

Starting with £/(3>) = f/o4> 7 , 7 _1 = wwo + § (we have chosen the positive square 

root by convention) and choosing the parameters Xu = 7 we obtain the potential 


U = 


l 2 Up <? 2 

4 — 67 2 


In order for a to be real we have chosen ojq > 0 which implies 0 < 7 < 3 . 

Introducing the parametrization [76] 

a = cr -1 exp(w + v), 
a = exp(7l(u — v)), 


(2.7a) 

(2.7b) 


where A is a constant to be specified, the Friedman equation (2.3c) for p^ = 0 becomes 

( 2 A 2 - 37 2 ) u 2 -2 ( 2 A 2 + 3q 2 ) iiv + ( 2 A 2 - 3^) v 2 + 7 2 U 0 = 0 . ( 2 . 8 ) 


Choosing the constant A = ^ 7 , (2.8) transforms to (see, e.g., similar equations (28) in [76] 

and (2.24) in [78]): 


uv = 


Uo 

12 ' 


Substituting the expressions 


Uqu . Uq 
V ~~12 u 2 ’ V ~ 12 u’ 


(2.9) 


( 2 . 10 ) 


the Raychaudhuri equation (2.3d) becomes 

3q 2 (l2u 2 + Uo) 2 + 4\/67 (l2u 2 + Uo) (—3 u — 12u 2 + Uo) 
+6 (l / 0 - 12 u 2 ) (-4 u - 12u 2 + U 0 ) = 0 , 

and the equation of motion for the scalar field (2.3a), now reduces to 

7 3 (I 2 u 2 + U 0 ) e VE ^ u ~ v '> (-4 y/Eu + 12 (7 - y/6) u 2 + (7 + y/6) U 0 ) 


= 0 . 


( 2 . 11 ) 


( 2 . 12 ) 


96 ( 37 2 - 2) u 2 

Since 7 is nonzero, it follows that both equations are simultaneously satisfied if and only if 

12 (7 — v/ 6 ) u 2 + (7 + \/ 6 ) Uo 


u = 


4^6 


(2.13) 
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The equation (2.13) admits the general solution : 


u(t) = < 


C 2 


21n(cosh(^^ (24 ^ +t >) 
y/6y —6 


, 7 2 <6 


2 In cos 


C 2 - 


V 7 2 -6 A /Uo(24ci+t) N \ 
2\/2 


-v/67—6 


—, 7 2 > 6 . 


(2.14a) 


Substituting the result for u on the equation (2.9), and integrating out the resulting equation 
we obtain 

2 ln(sinh ' 

C3+ \7 &t +6 " 


v(t) = 


7 2 < 6 


2 In ( sin ( Vr^ym^i+ V)) 


(2.14b) 


C 3 + 




272 


-v/67+6 


)) 


, 7 2 > 6 


respectively, where ci, C 2 and C 3 are integration constants 

rg 

3 


by hyperbolic functions: 


Since we have chosen 0 < 7 < given that ujq > 0, hereafter we use the branch given 


u(t ) = c 2 — 


2In fcosh ( V^VWaic 1+ t) \ 


2 In 


v(t) = c 3 + 


-v /67 — 6 

(sinh( ^^+‘> )) 


7/67 + 6 

Substituting (2.15) in (2.7) we obtain the solutions: 


cr(t) = eV 27 ^ t2 t3 ^ sinh ^ 7+6 (A(t)) cosh v^y-6 (A(t)), 

/IT ( \ 1 1 7 / 67+2 7 / 67—2 

a(f) = ev 2 ” C3 C2 ' C2 C3 sinh 7 ^ 7+6 (A(r)) cosh^ 7-6 (A(r)), 
-v/t7ocsch(2A(t)) ( 27/67 — 3 ( 7 s — 2) cosh(2A(i))) 


7/67 


\/b7 


#(*) = 


37/2^6 - 7 2 


where 


A(t) = 


7/6 - 7 2 t/Cq ( 24 ci + t) 

2t/2 


2.1.1 Including a massless scalar field 

The equation of motion (2.3b) for a massless scalar field is given by 


(2.15a) 

(2.15b) 


(2.16a) 

(2.16b) 

(2.16c) 

(2.17) 


a 

1 + 3—/> = 0, 
a 


(2.18) 


and it admits the solution 0 
parametrization (2.7) with A 



where is an integration constant. Combining the 
the Raychaudhuri (2.3d) and the equation of motion 
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for the scalar field (2.3a) we obtain 


(\/67 + 2) ( 37 s — 2) fi 2 exp ( 2 \/ 67 (u — v) — Gu — Gv) 

“= 37 

+ {'ll 1 ~ 3 ) + +3 ) ^ ~ (2 ' 19a) 

(•v /67 — 2 ) ( 37 s — 2 ) ,u 2 exp ( 2\/67 (u — u) — 6 tt — 6 u) 

+ ( 3 “ 3 \/l 7 )' i ” + (-\/l 7_3 )' >2 + 5^' (2 ' 19b) 

and the Friedmann equation (2.3c), assuming that V(4>) = 0, i.e., the massless case, now 
becomes 


uv = G(u, v ), 


G(u, v) = — 
v ; 12 


1 / (2 — 37 s ) /i 2 exp (2\/67 (u — v) — Gu — Gv) 


7 


+ Uq 


Combining the above equations we obtain 

3 (\/67 + 2) ( 37 s — 2) ^r 2 exp (2 \/67 (it — v) — Gu — Gv) 


247 u + 127 ^a/67 — 6^ u 2 + 7 67 + 6^ ? 7 o = 0, 


( 2 . 20 a) 

( 2 . 20 b) 


( 2 . 21 ) 


which reduces to (2.13) for /1 = 0. 

Now we want to simplify further the equations, choosing a new time parameter r such 

that 

u'v'f 2 = G(u,v), ( 2 . 22 ) 

where the comma denotes derivative with respect the new time r. Thus, choosing t = 
\JG(u, v) we obtain 

uv = 1. (2.23) 

The second derivatives with respect to t are given by 


u = t 2 u" + tv! , (2.24a) 

v = t 2 v" + tv', (2.24b) 


where 


r = 


/ (2—37 2 )/i 2 e 2v/ ®'7“—6( u +”) 

v -r 1 - + t/ ° 


r = — 


2^3 

(3 7 2 - 2 ) /r 2 ((V 67 - 3)- (V 67 + 3) v') e 2V&r(«-w)-6(«+v) 

12 ^ 


(2.25a) 

(2.25b) 
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Finally, the equation (2.21) transforms to 


u" (2 (3 7 2 - 2) M 2 e 2 ^“ - 2 7 2 [/ 0 e 6 “ +2 (^+ 3 ) 1 ') 

+ 7 (v^ (3 7 2 - 2) n 2 e 2 ^ u (u' 2 + l) + jU 0 (yfth + (y/Ey - 6 ) u' 2 + 6 ) e 6u+2 ( V ® 7+3 ) ,; ) = 0. 

(2.26) 


As in the previous case, since we are interested in the range of parameters 0 < 7 < y 
we omit the solutions of (2.26) given in terms of trigonometric functions and we use hyperbolic 
ones instead. Thus, for the unmodified case /i = 0we recover the exact solution (2.15) 


. 2 In (cosh (A (r))) 

“ (t) = - 

v(t) = c 3 + (siiih (A(r))) 


(2.27a) 

(2.27b) 


where A(r) = — 7 2 ( 2 ci +r), defined for 7 2 < 6 . 

We can use the solution (2.29) for constructing an approximated solution for the system 
when /i / 0 is a small parameter, i.e., assuming that the scalar field (f> is massless and has a 
small velocity cj). 

Assuming /r 7 ^ 0, we define 


u(t) =U(t) + fid 2 {r) + 0(f! 2 ), 

(2.28a) 

v(r) = V(t) + nd 3 (r) + 0(/i 2 ), 

(2.28b) 

where U, V are the seed solutions when fj, = 0 given by 


TTf n 2 In (cosh (A(r))) 

U(T) = C2 Ve-r-6 ’ 

(2.29a) 

V — c 3 + 2 lri (sinh (A(t))) 

[) 3+ ^ 7 + 6 ’ 

(2.29b) 


and d 2 ,d 3 are functions to be specified. Substituting in (2.28) and in (2.26), expanding in 
Taylor’s series with respect to the parameter /r near fi = 0, we obtain respectively: 


£11 + fi£ 12 + 0(/i 2 ) = 0, (2.30) 

£21 + 22 + 0(/i 2 ) = 0 , (2.31) 


where the equations £n = 0 , £12 = 0 , £21 = 0 , £22 = 0 must be satisfied. 


£11 = 0 => 

■ U'(t)V'(t) -1 = 0, 

(2.32a) 

£12 = 0 => 

■ U r (t )d 3 (t) + V' {r)d 2 {t ) = 0, 

(2.32b) 

£21 = 0 =4. 

■ VE 7 - 2 U"(t) + (y/E 7 - e) C/'(r ) 2 + 6 = 0 , 

(2.32c) 

£22 = 0 => 

• — 2d 2 "{T) + 2 ^-v /67 — 6 ^ U'(T)d 2 '{T) 


+ (2 (ve 1 

+ 3^ d 3 (r) + 6 d 2 (r) S j — 2U"(t) + ^v/fry - 6 ^ U\t) 2 + 6 ^ 

= 0. (2.32d) 
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Substituting (2.29) in the above equations makes the equations £\ \ = 0 and £21 = 0 trivially 
satisfied, and the equations £\2 = 0 and £22 = 0 simplify now to 


<h"(r) = 
dzij) = 


—A/fh/6 — 7 2 d 2 / ( t ) tanh (A(r)), 
(a/6 - 7 ) (fa 1 ( t ) coth 2 (A(r)) 

7 + Vg 


where A(r) 


-^\/6 — 7 2 ( 2 ci + r), with solutions 


d 2 (r) 

d3(r) 


3-1 1 


tanh (A(r)) 


\/ 6 “ 7 2 


+ /2j 


6 — 7/1 coth (A(r)) 


(7 + \/ 6 ) 


3/2 


+ /3; 


(2.33a) 

(2.33b) 


(2.34a) 

(2.34b) 


where /i,/ 2,/3 are integration constants. Henceforth, we obtain the first order (in the pa¬ 
rameter ^) solution 


. , 2 In (cosh (A(t))) 

“ M = C2 - Ur -a ^ 


, . 2 In (sinh (A(r))) 

w(r) = c 3 +- V 7 = v + n 


3/1 tanh (A(t)) 


a/6-7 2 


+ /2 


+ o (^ 2 


1/67 + 6 


6 - 7/1 coth (A(r)) 


(7 + a/ 6 ) 


3/2 


+ /3 


, (2.35a) 

+ 0(/r 2 ). (2.35b) 


The relative errors in the approximation of (2.35) by (2.29) are: 

\/f/i tanh(A(r)) 


E r (u(r )) : = 


E r (v(r)) := 


u ( t ) — C/(r) 

«(r) 


u(r) — V (r) 
u(r) 






+ /2 


C2 - 


2 ln(cosh(A(r))) 
V67—6 


+ o (/, 2 ) , 




3 v A- 7/1 coth(A(r)) 

- + 73 


(7+7/6) 


3/2 


2 ln(sinh( A(r))) 

7/67+6 ^ 63 


+ O (/i 2 ) 


(2.36a) 


(2.36b) 


Taking the limit r -+ +00 it follows that the above relative errors tend to zero. Thus, the 
linear terms in [i in the equation (2.35) can be made a small percent of the contribution of the 
zeroth-solutions (2.29) taking r large enough. Henceforth, this shows that the behavior of 
the solutions for the induced gravity model does not change abruptly when a small massless 
scalar field, 0 , is added to the setup. 
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Finally, going back to the original variables, we obtain the solutions: 


7/67 


7/67 


ct(t) = e* 63 ^ sinh T6 7 +6 (A(r)) cosh 76 7 -6 (A (r))x 

1 ^ ^ j. ^ 4 /icsch( 2 A(r)) (7/6- 7 cosh( 2 A(r))) 

1 + 2 7 M V6(/ 2 - /s)- i ^ - 


7/7/6 - 7 (7 + 7/6)' 


+ 0(/i 2 


( 2 . 37 a) 


o(r) = 


7/67+2 7/67—2 

sinh 767+6 (A(r)) cosh^T- 6 (A(r)) 


A' 


1/3 


+ 


+ 


fiH ( 2 7/67 — 3 (7 s — 2) cosh( 2 A(r))) sinh v^7+e (A(r)) cosh^T- 6 (A(r)) 

37 / 7 / 6-7 (7 + \/ 6) 3/2 ^i 1/3 

/2/i (7 (7 (7/67 + 10) + 2 7/6) — 12) sinh( 2 A(r)) sinh ^67+6 (A(r)) cosh(A(r)) 

4 (7 + VS ) 2 i^i 1/3 

fsH (7 (7 (7/67 + 14 ) + 10 \/ 6 ) + 12 ) sinh( 2 A(r)) sinh Ve-t+a (A(r)) cosh 767-6 (A(r)) 9 

-h 0(/i ), 


4 ( 7 + \/6) 2 A' 1 1/3 


( 2 . 37 b) 


A(r) = 


fa'(r) _ t/G(«(t),?;(t)) a'(r) 


+ 


+ 


+ 


+ 


o(r) a(r) 

Y / A(r)csch(2A(r)) ( 47/37 — 3\/2 ( 7 s — 2) cosh(2A(r))) 

67 / 6 — 7 2 

27 / 1 /i coth(2A(r))csch(2A(r)) ((30 — 97 s ) A(r) + ( 77 s — 18) C/o) 
7/3 (7 + 7 / 6 ) ( 7 2 - 6 ) \J F(t) 

\/2/i/icsch 2 (2A(r)) (( 27 4 — 97 2 + 18) F(t) + (- 7 41 + 7 2 — 6 ) C/ 0 ) 
(7 + 7 / 6 ) ( 7 2 - 6 ) yjF{r) 

7/2 (7 4 — 57 s + 6) /i//cosh( 4 A(r))csch 2 ( 2 A(r))(t/o — F(t)) 

(7 + t/6) ( 7 2 - 6) 7 /F(t) 

a 2 /icsch( 2 A(r))(A(r) — Uq) ( 7/6 ( 7 2 — 2) cosh(2A(r)) — 47 ) 

27t/A(t) 

cc 3 /i(A(r) — C/ 0 ) ( 7/6 ( 7 s — 2) coth(2A(r)) — 47 Csch( 2 A(r))) 


= 4> 0 + 


/ 

Jtc 




4 7t/A(t) 
-dr’ 


( 2 . 37 c) 


= />o + 




'T 0 


£ — + 


f 

J TQ 


j 27 / 3 /i A/ sinh 767+6 3 (A(r , ))cosh T6 7 -6 3 (A(t / )) 

\ 7W) 

dr' f T 27/3 , . 

—= to + , dr 

7 /G(u(t'),v{t')) Jtq 7 /F(t') 


dr' + 0 (/i 2 ), ( 2 . 37 d) 


+ 


+ 


r 

J TO 

r 

Jtq 


2 /i/icsch(A(r / ))sech(A(r , ))([/o — F(t')) ( 2 7/3 (7 s — 3 ) cosh( 2 A(r / )) — 37/27) f 


7/7/6 — 7 (7 + 7/6)A(r') 3 / 2 


6 ( 7/3 — 7/27) / 2 /i(A(r / ) — 1 /p) 
A(t') 3 / 2 


dr'+ f T 

J TQ 


6 (7/27 + t/ 3 ) / 3 /i(A(r') - C 7 o) 


A(r) 3 / 2 


dr 

+ 0 (/i 2 ), 

( 2 . 37 e) 
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which generalize solutions (2.16). 

In (2.37) we have introduced the expressions 


K l — e 3 \f^'y( c 2-C3)-3(c2+C3) 



\ 2 

C] j e 2v/67(c2-c 3 )-6c2-6c 3 


OL 2 = 


\/\/6 - 77 (\/ 27 3 + 3 \/ 37 2 - 6\/3) / 2 
(7 + a/ 6) 37,2 ( 7 2 - 6 ) 


a 3 = — 


2 V / \/6-77 (\/27 3 + 5V37 2 + 12x/2 7 + 6>/3) /s 

/7. , ,/7t 3 / 2 ^. 2 ^ 


(7 + \/6) 7 (7 2 — 6) 


and 



We note that by fine-tuning the parameter values we may get some integrable cosmological 
models in the case of a massless and slowly moving scalar field. This makes our choice of 
potential very interesting. This issue, and the discussion of massive scalar field, q i, deserves 


further investigation and it is left to future projects. However, it is worth noticing that the 


main focus of this research is not to find analytical solutions but to study the asymptotic 


behavior of the solutions space without using fine-tuning of the parameters and the initial 
conditions. Dynamical systems theory is a powerful tool for doing this research. 


2.2 Dynamical system analysis 

In order to study the cosmological behavior in a general way, independently of the initial 
conditions and the specific universe evolution, we apply the dynamical systems method, which 
allows to extract the global features of a cosmological scenario [80-95] . In this procedure, one 
first transforms the involved cosmological equations into an autonomous system and then one 
extracts its critical points. Hence, taking linear perturbations around these critical points, 
and expressing the perturbations in terms of a perturbation matrix, allows to determine 
the type and stability of each critical point by examining the eigenvalues of this matrix. 
In the case of non-hyperbolic critical points one should use the center manifold theorem 
[82, 83, 86 , 88-91, 93, 96, 97], 

2.2.1 Finite analysis 

Let us define the following dynamical variables 


e = n/T, 



(2.38) 


The Friedmann equation (2.3c) leads to 



(2.39) 


It is defined the auxiliary variable [26]: 



(2.40) 
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that it is interpreted as the dimensionless kinetic energy density of the JBD held 4>, and it is 
not necessarily positive. This implies that x 2 + z 2 might be greater than the unity. This is 
due to is a non-conventional scalar held that does not follow the standard energy conditions 
for a scalar held in GR. Besides, physical conditions > 0, U$ > 0 and H > 0 implies z > 0. 

The evolution equations for the variables (2.38) are given by: 


x = 3 - 


9y 2 


x 3 +\ye - 




+ 


+ 3 (y 2 - l) x + 


(3y 4 — 5y 2 + 2) xy 2 
4y 2 


1 


+ 2 ( 6 7 2 -5 )xy + -7 xz 2 {\u - 2 y), 


y = 3y 2 (2 - 3x 2 ) + y ^9y 2 + ^3 - x 2 + ^ yz 2 (Xjj - 2 y) - 3^j + 


+ 7 (zr + ^ - 5 ) r + 


97 2 


«' = * ^3 7 2 + ^3 - ® 2 ^ + ^ (V/ 


- 4 y 2 + 3'yz 2 {\ u - 27 ), 


2 + — - 5 ) y 2 z+ 


V 


e = 


yz (- 67 s + 37 + Xu) 3 3 

-^-+ ip z (hi - 27 ), 

yj_ 

2 ’ 


(2.41a) 


(2.41b) 


(2.41c) 

(2.41d) 


where the comma denotes derivative with respect the conformal time r = In a. 

From the equations (2.41) it follows that the signs of e and z are invariant in time. This 
means that all the solutions with e(ro) = 0 (respectively e(ro) < 0 or e(ro) > 0 ) at an initial 
time ro, will satisfy e(r) = 0 (respectively e(r) < 0 or e(r) > 0) at any time r. Given that 
the sign of e is invariant for the how, we can safely analyze the stability of the hxed points 
with e = 0 in a neighborhood of it, containing both e > 0 and e < 0 points, but since the 
boundary e = 0 cannot be crossed, we restrict ourselves to the region of physical interest. 
The same reasoning applies to z. The “recipe” to deal with this kind of dynamical systems 
was given in the seminal work [81]. Henceforth, we can focus our analysis on solutions with 
e > 0 and z > 0, and the equations (2.41) define a flow on the phase space: 


^ := [(x,y,z,e) € M 4 : x 2 + y (^-y - l) + z 2 < 1 , e > 0, z > oj . 
Furthermore, the cosmological parameters are given by: 

_ pi _ 2 y 2 e 2 + 3q 2 (4 - e 2 {{y + 2 ) 2 - 4z 2 )) 


Hi = 


Qo — 


3 H 2 
P2 

3 H 2 


12 7 2 

(3 7 2 ((y + 2 ) 2 - 4z 2 ) - 2 y 2 ) 


8q 2 


H 


9q 


2 \ 1 

,2 1 


q = -1 - --7 = -1 + 3Y + 3 - -2- )x z + - [3Y + -r - 5 + 


H 2 


+ (37 2 - 2 ) y + ~'yz 2 {X u - 27 ), 


7 


(2.42) 


(2.43a) 

(2.43b) 


(2.43c) 
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2 H 


Wtot = -1 - = -! + 2 7 2 + (2 - 37 2 ) 3:2 + p ( 3 7 2 + 72 “ 5 ) y 2+ 

+ ( 2 7 2 - 7 ) y + 7 z 2 {\u - 27 ). 


(2.43d) 


2.2.1.1 Fixed points and stability in the Jordan frame. The critical points of the 
system in the Jordan frame (2.41) are: 

1. Ji ,2 • (x = 0, y = ± , z = 0, e = 0^ . They always exist. The eigenvalues are 


37 


37T\/6’ 


3 7 6(V6t2 7 ) 3(V6tAc/) \ 

3 7=pV6’ \/6=F 3 7 ’ \/6=F 3 7 )' 


J\ is a sink for y/| < 7 < y |, Xjj < 77 
saddle otherwise. J 2 is always a saddle. 


source for A u < x/ 6 , 0 < 7 < 71 or a 


2. J 3 : ^0, ,0; 0 ) ■ It exists for —\J | < 7 < — 1, or— 1 < 7 < 1, or 1 < 7 < The 

eigenvalues are 

f 7 2 0 3 - 27 s 7(A[/ - 27) ] 

1 7 2 -l’ ’ 7 2 -l 5 7 2 — 1 J ' 

It is always a saddle. 

3. J 4 : ^0, 0, ^==, 0^ . It exists for 7 > 0, A{/ < 0. The eigenvalues are 


jo, -3, i (-3 - V / 24 7 A c/ + 9) , i (-3 + V247A u + o) j 


It is nonhyperbolic with a 3D stable manifold for A u < 0,7 > 0. Thus, it has a large 
probability to attract the universe at late times. The full stability analysis requires the 
application of the center manifold theorem (the analysis is done in subsection 2 . 2 . 1 . 2 ). 


4 - ^ 5,6 



always exist. The eigenvalues are 


j -1, 2, i (-v 7 ! - 24 7 2 - 1) , \ ( s/\ - 24 7 2 - l) } . 


Thus, they are always saddle. 

5. J 7 : ^0, - $u -2 > ‘ 11 existsfor -| < < 0,-^6 < Xu < 72V2W0 + 3 

or —| < wo < 0,7272^+3 < A;/ < 76 or wo = 0,-76 < Xu < 76 or wo > 
0, —76 < Xu < 76. The eigenvalues are 

f 7 Aq 6 ~ A {/(7 + At/) 2A[/(27 — Xu) 6 — A u 2 \ 

1 2 - jXu ’ 7A[/ — 2 ’ jXu - 2 ’ 7A u — 2 J 

It is always a saddle. 
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6- ./s,9 : ( ± 


_3 ^ ^ _5 

2 < Wo < g. 


Vh 7 h+*u )~6 67 ^V7+Aq n \ Um , ov5of 

V(A7-2 7 ) 2 ' A ^’ V(Ao-2 7 ) 2 ’ V • 67 6 St ° 


w 0 


y /(\ u - 2 fy )* ' lu - VV ’ y / i A„ „ x 

5 . / 12 ^o+V^o+ 37+19 <\u<VQ, or 

= | (a/ 102 - \/6) < A;j < \/6, or < w 0 < - 5 , ^ 

°r "I < w 0 < -5,75^5 < A[/ < V6, or w 0 


1 


2\/2 

V 2 / 0+3 

W 0 > - 


| 12 ^°+V24ajo+37+ i9' < ^ < y^. xhe eigenvalues are 


12cjQ+V24a;o+37+19 ^ ^ 

2cjo+3 — U 

= — 5,2 < A (7 < \/6 or 


f 37 _Mi_M 2 _ M-3 ) 

\ A[/ - 27’ 2 ' y(Xu - 27) 3 ’ 27 (A ;7 - 27)2 ’ 27(A[/ - 27) 3 J ’ 

where mi,M 2 ,M 3 are the roots of P(/r) = -1447 4 (7 + Xu)(^u(l + Ac/) - 6)(2^ - Xu) 7 + 
127 3 m(A; 7 - 27) 4 (( 7 2 - 12) A[/ + 2^Xu 2 + 67 + A[/ 3 ) - 67 /i 2 (57 -2Xu)(Xu -2^) 2 + m 3 • 
Their stability should be analyzed numerically. 

The existence and stability conditions of the above critical points are displayed in Ta¬ 
ble 1, whereas the cosmological parameters (2.43) evaluated at the critical points, and the 
description of these critical points are given in Table 2. 


- 14 - 


























Label 

x, y, z, e 

Existence 

Stability 

J l,2 

(° ■ ± vg 3 ,> 0 ' 0 ) 

always 

Ji is a sink for — | < u>o < 0, A[/ < VG 

Ji is a source for Xu < VG, ojq > 0 

J\ is a saddle otherwise 

J2 is always a saddle 

Js 

(O.^r, 0,0) 

CJo ^ 2 ’ ^ 6 ^ ^0 ^ 2 

saddle 

Ji 

(°’ 0 ' J^’°) 

Wo > — A(7 < 0 

sink for cuo > — |, Ay / 0) A u < 0. 

^5,6 

( ±V /S,0,0,±£) 

Xu > 0 for J5 and Xu < 0 for 

saddle 

•h 

/n / 2 / 6-Ar/ 2 n\ 

W’ 7 Ac/- 2 > V 3 V ( 7 A C7 - 2 )^ U / 

—\/6 < A u < \/6,7 > 0 , 7 Ac/ / 2 

saddle 

Js ,9 [ 

' 1 y/^uil+^u )—'6 67 V 7 V 7 +AU n ' 

\n^ \., <r * fci - 6—A r/ \- L 0^, r,r 

numerical inspection 

v V( A ^- 2 7 ) 2 ' Ac/ ~^ 7 ' V(Au- 27 ) 2 ' "y 

1 u <. Af/ \ V 0,7 7i ^ ? A// 7= Z7 or 

7 > 0, A[/ > \/6, Af/ 1 27 


Table 1. The existence and stability conditions of the critical points of (2.41). We use the definition 7 = (wq + 3/2) 1 ' 2 . 

























Label 

n. 

fl2 

Df 

q 

^tot 

Description 

Hit) 

•h,2 

i 

0 

1 

2 V 6 

7±\/6 

./i 2 are dominated by f 2 ^. 

Ji is accelerating for — | < wo < 0. 

J 2 is always decelerating. 

Ho 

\/6=f3 7 

\/6t37 

1+3 ^(^)( 4 

h 

i 

0 

7 2 (3 7 2 -4) 

l-2 7 2 

3-57 2 


Ha 

3(1—T 2 )^ 

1 — 7 2 

3(1-7D 

UCIWCC11 3 -L'D-*- U/y -f— 

Dominated by for u>q = — |, 

Accelerating for wo > 5 or — | < wo < — \ ■ 


Ja 

i 

0 

0 

-1 

-1 

Intermediate accelerated a(t) ~ e“ ltP1 , ai > 0,0 < pi < 1. 

~ Q!iPit Pl_1 . 

•hfi 

1 4 

1 — T 2 “ 
A V 

4 

A v 

0 

1 

1 

3 

Radiation dominated. 

Ho 

l-\-2Ho(t—to)' 

■h 

l 

0 

■V/(-37(7-V/-4)-2Au) 

2 -A 2 , 

6 —Aj/( 7 + 2 A[/) 

Dominated by the energy density of $. 

Ho 

3( 7 Ae7-2 ) 2 

7 At/-2 

37A[/-6 

1 +" oA X l 7 2 Acr, (* <o)' 

^8.9 

l 

0 

3( 7 2 -2 7 Ar;+2) 

7-2A r j 

A[/ 

Scaling solution. 

Ho 

(Ac- 2 7 ) 2 

27 -A u 

27 - Air 

1 + 3 ff o( 1 ^ r )(t-to)' 


Table 2. Description of the cosmological parameters (2.43) of the critical points of (2.41). 




































2.2.1.2 Center manifold analysis for the intermediate accelerated solution J 4 . 

In order to investigate the stability of the center manifold for J 4 we introduce the new 
variables 


u = e, 

XuXye 

th = —=-h x, 

\/6(27 - A u) 

V7 (12 7 2 ~ V24 7 A u + 9-3) ^ y ( 6 7 3 - 3 7 - A v ) 

V12 7 — 6A{/v / 87A;7 + 3 y / 6 v / 7\/27 — + 3 

z (-12 7 2 + V24 7 A u + 9 + 3) 

+ 2V24 7 A u + 9 

V7 ( 12 7 2 + V24 7 A u + 9 -3) y (- 67 3 + 37 + Xu) 

Vo —- — - —|—-1— 

\A2 7 — 6\uV8l^u + 3 vM/tMt — Xu V ^7 Xu + 3 

z (12 7 2 + V247AC/ + 9-3) 

+ 2V247A{/ + 9 


(2.44a) 

(2.44b) 


(2.44c) 


(2.44d) 


to obtain 


f u '\ 

[ 

v'l 

= 

v'2 


\ V 'J 

V 


0 -3 
0 0 


0 

0 

\ 

/ u \ 


(f\ 

0 

0 


Vi 


9 1 

\/2AyXu + 9) 

0 


V2 

~r 

92 

0 

-4 (3 + ^/2A^Xu 

+ 9)/ 

\v 3 ) 


\W 


(2.45) 


where {f,gi,g 2 , 93 ) T is a vector of higher order terms. 

Since the center subspace of the origin is tangent to the e-axis, it follows that the center 
manifold of the origin is given locally by the graph 


{(u,ui,u 2 ,u 3 ) : vi = hi(u),v 2 = h 2 (u),v 3 = h 3 (u), 

Mo) = Mo) = Mo) = o, 

Mo) = Mo) = Mo) = 0 , M < < 5 }, (2.46) 


where 5 is a positive small enough constant. The functions hi,i = 1,2,3 satisfy a set of 
quasilinear ordinary differential equations which can be expressed symbolically as 


Ki u ) u ' ~ v i 


Vi=hi(u) 


= 0, z = 1, 2,3, 


(2.47) 


where one must substitute v!, v[, v' 2 , v' 3 through (2.44) and use the replacement v\ —>• h\(u),V 2 —>• 
h, 2 {u) and v 3 —>• h 3 (u). 

Setting 


h\(u) = anu 2 + ai 2 U 3 + O (u) 4 , (2.48a) 

h 2 {u) = a, 2 iu 2 + C 122 U 3 + O (u ) 4 , (2.48b) 

h 3 {u) = a 3 iu 2 + a 32 u 3 + O (u ) 4 , (2.48c) 
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in (2.47), equating to zero all the coefficients of equal powers of u, and solving for the a*/s 
we get up to fourth order: 


an — 0, ai2 — 


Ac/Ay( 7 Ac/ - 3) 


021 = 


«31 = 


3>/6(At/ - 2 7 )3 ’ 

(7 (3 7 v / 8 7 A Lf + 3 + 5^7 - 2 V3Xu) - 2 (/8 7 A v + 3 + V3)) 
21/2(27 - Ay) 5 / 2 (8^ 7 Ac/ - 3/87AC/ + 3 + 3^3) 
y^fX'ljXy (7 (-37V87AC/ + 3 + 51/37 - 2i/3Ac/) + 2/8 7 Ac/ + 3 - 2^3) 


, 022 — 0, 


2^(27 - Ac/) 5 / 2 (8 a/3 7 Ac/ + 3 (V8 7 Ac/ + 3 + v^)) 
Henceforth, the dynamics on the center manifold is given by 


u = 


l^Xl 


u 


+ O {uf 


2(A u - 2 7 ) 2 

Neglecting the hfth-order terms and integrating we hnd that 

\J-{Xu - 2 7 ) 2 


a(r) = ± 


^2d{Xu - 2 7 ) 2 + 'yXuXyT 


(2.49a) 

(2.49b) 


, 032 — 0. (2.49c) 


(2.50) 


(2.51) 


where ci is an integration constant that must be negative in order for u to be real-valued. 
Thus, for 7 > 0, Xu ^ {0,2 7 },Ay / 0, it follows that the origin, and then J 4 , is stable 
provided Xu < 0 . 


2.2.1.3 Special case: A u = 0, [/($) oc <h 2 . We introduce the new variables 


u 1 = e, 

o 2 = ^(2 7 2 - 1) fo-2z + 2), 


ui = x, 

«2 = | -7 2 (y- 22; + 2 )- 


( u>1 \ 


/0 0 0 0 \ 


/wr\ 

u' 2 


0 0 0 0 


U2 

v[ 


0 0—3 0 


Vi 

\ v U 


CO 

1 

O 

O 

O 


\V 2 ) 


+ higher order terms. 


(2.52a) 

(2.52b) 

(2.52c) 

(2.52d) 


(2.53) 


Since the center subspace of the origin is tangent to the plane U 1 -U 2 , it follows that the 
center manifold of the origin is given locally by the graph 


U2 ,vi,v 2) : vi = hi(ui,u 2 ),v 2 = h 2 (u l ,u 2 ), 
hi(0,0) = h 2 (0,0) = 0, Dh(0, 0) = 0, u\ + u\ < $}, 


(2.54) 


where Dh is the matrix of derivatives and 5 is a positive small enough constant. The func¬ 
tions hi, h >2 satisfy a set of quasilinear partial differential equations which can be expressed 
symbolically as 


~ dhj(u) 
. du\ 


u x + 


dhj{u) 

du 2 


«2 - V i 


Vi=hi(u) 


= 0,z = 1,2, 


(2.55) 
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where one must substitute u\. u ' 2 , v \, v ' 2 , through (2.53) and use the replacement vi —>• 
hi(ui,u 2 ),v 2 ->• h 2 (u 1 ,u 2 ). 

Setting 


hi = auu\ + a\ 2 uiu 2 + a 22 ul + 0(3), (2.56a) 

h 2 = buu\ + bi 2 uiu 2 + b 22 ul + 0(3), (2.56b) 


and plugging back in (2.55), equating to zero all the coefficients of equal powers of u\ and 
u 2 , and solving for the a^-’s and b^ ’s we get up to third order 


«11 «12 cl 22 
b n b 12 b 22 


( 


0 


2 7 1 2 3 -l 


0 


o o t+++ 


3(l-2 7 2 )^ 

Thus, the dynamics on the center manifold is dictated by 

27 2 «irt 2 


ui = 


u 2 = 


1 - 27 2 
87 2 u 2 


+ 0(3), 


1 - 2Y 


+ 0(3), 


(2.57) 


(2.58a) 

(2.58b) 


where 0(3) denotes error terms of third order in the vector norm. Neglecting the error terms 
and integrating out the system (2.58) we obtain 


ui = c 2 v/ 7 2 ( 8 r — 2 ci) + ci, u 2 


2 7 2 - 1 

7 2 ( 8 r — 2 ci) + ci ’ 


(2.59) 


Finally, it follows that for A u = 0, the origin, and then the point J 4 behaves as a saddle since 
the orbits departs from the origin along the e-direction as the time goes forward. 

2.2.1.4 Features of the critical points of the system (2.41). Let us summarize the 
features of the critical points of the system (2.41) found in subsection 2.2.1.1: 

1 . J 72 is dominated by the kinetic term of <L, that is 17^? = 1 , and the quintessence 

held has Vt 2 = 0. J\ is a sink for — | < ujq < 0. This range for wo is several orders 

of magnitude lower than the bound ujq > 4 x 10 4 * * imposed by the Solar System tests 
[22, 23], the bounds estimated on the basis of cosmological arguments ujq > 120 [24] 
and 10 < wo < 10 7 [25]. Therefore, in section (2.3) we will discuss on the asymptotics 
of Ji when it is a source, i.e., for A jj < V6,ujq > 0. J 2 is always a saddle. 

2 . J 3 is a scaling solution between 17^ and for ujq / ~~\[\- this case, —1 

9 (rJy 5)2 an d —>• 0. It is dominated by f7^ when ojq = — y^f, and accelerating when 

wo > \ or — | < ujq < — 

3. J4 represents an accelerating solution with rctot = — 1. At the end of section 2.3 we 

will discuss on the corresponding asymptotics giving new arguments supporting the 

statements in [26] against [19-21]. Indeed, J4 represents an intermediate solution at 

late time, i.e., a(t) ~ e aitP1 as t —>• 00 where 07 > 0, and p\ := — Al 2 _ 2 37 , 0 < pi < 1 

provided 7 > 0, Xjj < 0. Additionally, rc to t —>• — 1 as t — > 00 for 7 > 0, Ay / 0, Xjj < 0. 
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4. represent a radiation-dominated solution (u>tot = 7 ) which are always saddle. 

5. J 7 denotes a solution dominated by the energy density of the scalar field <f>. It is always 
a saddle. 

6 . Js : 9 represent scaling solutions where the contributions of and </> to the total energy 
density are of the same order of magnitude. 


2.2.1.5 Hubble parameter of solutions near the fixed points. Now, we obtain 
first order approximation for the Hubble parameter of solutions near the fixed points of the 
system (2.41), the results are presented in the last column of Table 1. By definition we have 
H = —(1 + q)H 2 , with q defined by (2.43c). This expression is valid in the whole phase space 
and not just at the fixed points. In the case when q 7 ^ —1, by continuity we can approximate 
the value of q for a solution close to a given fixed point by the constant value q = q* 7 ^ — 1 , 
where the asterisk means evaluation at the fixed point. 

Defining the reference values Hq = H\t=t 0 , where Hq and to are finite numbers, we 
obtain 


H(t) 


Hp 

1 + H 0 (l + q*)(t - t 0 ) ’ 


(2.60) 


which provides a first order approximation for the Hubble parameter at a given fixed point. 

These arguments are valid for all the fixed points, with the exception of J4 since q = — 1 
for it. For the point J4 the above procedure can not be applied, instead we proceed as in 
Section 2.3 to obtain H ~ «iPit pi_1 , with 0 < p\ < 1. Additionally, for the fixed points 
which are saddle this approximation is valid whenever the stable manifold of the fixed point 
is approached 1 . 

From Table 1 we note that the fixed points J 4 , J 7 , Jg and Jg have e = 0 (i.e., <!> = 0) 
and a finite non-zero value of z oc 4> P /H, where p = -y ~'* u . This means that H and <3? p ( p 7 ^ 0) 
both tends to 00 or to 0, at the same rate, depending on the sign of p. That is, for 7 > A u 
we have 4> p —>• 0, H —>• 0 and for 7 < Au we have 4> p —>• 00, H —>• 00. On the other hand, 
from Eq.(2.60) it follows that either H —>• 0 as t —>• 00, or H —>• 00 for a finite value of time 
t —> 1 which corresponds to a finite time singularity (see references [98, 99] and 

the more recent work [100] for the classification of finite-time future singularities). Thus, we 
can combine the above facts to determine what happens with H at J 7 , Js and Jg, at least 
asymptotically. 


1. For J7 we have (for solutions starting at the stable manifold of J7): 

(a) — \/6 < A[/ < 0 ,7 > 0, H —>• 0. 

(b) 0 < A u < \/2, 0 < 7 < A jj, H —> 00 (at finite time). 

(c) 0 < A u < y/2, X(j < 7 < H —y 0. 

(d) 0 < Ajj < \/ 2,7 > H —>• 0. 

(e) y/2 < Au < \/ 6 , 0 < 7 < H —> 00 (at finite time). 

(f) y/2 < Au < \/ 6 , < 7 < Au, H 00 (at finite time). 

(g) V2 < Au < x/ 6,7 > Au, H —>• 0. 

x For the solutions starting at the unstable manifold of the fixed point, the orbits depart from it, and the 
expression (2.60) might not be accurate. 
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2. For Jg^ we have (for solutions starting at the stable manifold of Js, 9 )- 


(a) 

0 < Xu < \/ 3 ,7 > 

,H^ 

0. 


(b) 

V 3 < Xu< 2 , 6 -^ <7 

< Ac/, 

H -> 

00 (at hnite time) 

(c) 

2 < Xjj < y/ 6 , < 7 

<r All 
' 2 ’ 

H -> 

00 (at finite time) 

(d) 

2 <Xu <y/Q, ^ < 7 < 

Xu, H —> oc 

(at hnite time). 

(e) 

< A u < VE, 7 > Xu, 

H -> 

0. 


(f) 

Xu > VG, 0<7<^p,i7 —»oo (at finite time). 

(g) 

A u > VG, 77- < 7 < Xu, 

H ^ 

00 (at 

hnite time). 

(h) 

Xu > VG, 7 > Xu, H —»• 

0. 




For J 4 (2.60) does not apply. However, since this point exists for 7 > 0 and A u < 0, 
and it is stable for 7 > 0, Xu < 0, then p > 0, thus, as <f> —>• 0, H —>• 0 as the stable point is 
approached in order to z be a non-zero constant at the fixed point, as verified in Sect. 2.3. 


2.2.2 Analysis at infinity 

As stated before is not necessarily positive, therefore, x 2 + z 2 might be greater than the 
unity, which implies that the phase space (2.42) is not compact. So, in order to obtain global 
results about the dynamics, we implement a compactification scheme. For this analysis, the 
following compact variables are useful: 

X = -,Y = ~, Z = -, r= ^l + x 2 + y 2 + z 2 , (2.61a) 

r t v 

and 

E = jl-, (2.61b) 

which satisfy the evolution equations 


dX _ EX V (X 2 - 1) [12 7 2 + 127 2 KY - 24 7 2 A 2 - 97 2 Y 2 - 2 V 2 - 24 7 2 Z 2 ] 
dT 4-\/67 2 

(E - 1 )KX [6 7 2 ( 7 2 (5A 2 + 4Z 2 - 2) - 4A 2 + 2) + (33 7 4 - 7 7 2 - 2) Y 2 ] 

+ Yf + 

(E - 1)A uXZ 2 [37 2 A - 6 7 2 Y + Y] ( 

27 + 

1 )XY [37 2 (8A 2 + Y 2 + 6 Z 2 - 2) + 5A 2 + 8Y 2 + 3Z 2 - 5] , (2.62a) 

dY _ EX V XY [37 2 (-4 KY + 8A 2 + 3 Y 2 + 8 Z 2 - 4) + 2Y 2 ] 
dT 4\/67 2 

(E - 1 )KY [6 7 2 ( 7 2 (5A 2 + 4Y 2 - 6) - 4A 2 + 2) + (33 7 4 - 7 7 2 - 2) Y 2 ] 

+ Yf + 


(E - 1)A uZ 2 [3 7 2 (Y(A - 2Y) + 2) + Y 2 ] 

27 + 

-\{E- 1) Y 2 (5A 2 + 8 Y 2 + 3Z 2 - 8 ) + 

+ 37 2 (2A 2 (4Y 2 - 5) + Y 4 + Y 2 ( 6 Z 2 - 5) - 8 Z 2 + 4) 


(2.62b) 
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dZ _ EX V XZ [3 7 2 (-4 KY + 8 Y 2 + 3 Y 2 + 8Z 2 - 4) + 2Y 2 ] 
dT 4\/6 7 2 

(E - 1 )KZ [-24 7 2 V 2 + 6 7 4 (5X 2 + 4Y 2 - 2) + (33 7 4 - 7 7 2 - 2) Y 2 ] 

+ 4^ + 

(E - l)\uZ [~3^ 2 Z 2 {K - 2Y) - YZ 2 + Y] 

+ 2 ^ + 

- ^(F - 1)YY [3 7 2 ( 8 Y 2 + Y 2 + 6 Y 2 - 2) + 5X 2 + 8 Y 2 + 3Y 2 - 3] , (2.62c) 

r] F 1 1 

— = - Y(E - lfE, (2.62d) 

where K = \/l — X 2 — Y 2 — Z 2 , and we have introduced the new time variable T given by 
dT = (1 — E)~ 1 K~ 1 dr. The physical region of the phase space becomes 

Too := | (X, Y, Z, E) : 2X 2 + Y 2 (y + l) - Y\J\ - X 2 - Y 2 - Z 2 + 2Y 2 < 1, 

Y > 0, Y 2 + Y 2 + Z 2 < 1,0 < E < l}. (2.63) 

The critical points of the system (2.41) at infinity are the critical points of (2.62) located on 

|(X,Y, Z,E) : 2X 2 +Y 2 (y + l) - Yy/l - X 2 - Y 2 - Y 2 + 2Y 2 < 1, 

Y > 0, Y 2 + Y 2 + Y 2 = 1 , 0 < E < lju 

| (X,Y,Z,E) : 2X 2 +Y 2 (■y + i) - Y^l - A ' 2 - Y 2 - Y 2 + 2Y 2 < 1, 

Y > 0,X 2 + Y 2 + Y 2 < 1,E = lj. (2.64) 

The critical points of the system (2.41) at infinity, located at the Poincare sphere X 2 + 
Y 2 + Y 2 = 1, are shown in Table 3. None of these points can be considered of physical 
relevance since their existence conditions are in stress with the observational bounds on 
ujq obtained in [22, 23, 25]. On the other hand, in order to describe the critical points 
at infinity, which satisfy E = 1 and X 2 + Y 2 + Y 2 < 1, we choose: the variables x,y,z 
defined in (2.38), the variable E defined in (2.61b) and a new time variable T defined by 
dT = (1 — E)~ 1 dr. The resulting system admits the set of non-hyperbolic equilibrium points 

2 at infinity z c = \J 1 — x 2 — + y c , E = 1, representing the boundary of the phase space 

T but with T —>• oo. 

2.3 Viability of the intermediate accelerated solution in the Jordan Frame 

The possible future attractors of our model, in the Jordan frame, are critical points J\ and 
J 4 . Ji is an attractor for — | < wo < 0, but this existence condition is at variance with 
the observational bounds reported in [22-25], thus, the discussion of it is omitted 3 . On the 
other hand, J 4 is a late-time attractor for 7 > 0, \y / 0, Xjj < 0, and the equation of state 
parameter becomes u; t ot = — 1 as the critical point is approached. Let us discuss more on 
the corresponding asymptotics. 

2 All the eigenvalues of the Jacobian matrix are zero. 

3 We restrict ourselves to positive values of BD parameter, wo > 0, despite of the observational results 

reported in [22-25] could be less restrictive in the cases of modifed JBD theories. 
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Label X, Y, Z, E 


Existence 


Ql,2 


( 0 , ± 1 , 0 , 0 ) 


— 4 < wo < 0 




2ll)q + 1 


2x/3 


2u 0 -ll> VH-2^0 


, 0,0 


o ± Wo < 9 


Q 5,1 


Q9, 


10 


2oin + l 


2±3 


2^0-11’ s/ll-2u 0 


, 0,0 


-| < Wo < -5 


^7,8 (o,y 7 * 8 yi-y 7 * 8 2 ,o) -I < wo < -i ^y 2 } < A,/< v / 4^±6 


± 


vF 

%/wo — 6 ’ y/6—U0 




, 0,1 


— 4 < wo < 0 


Q11, 


12 


± 


y/^0 V§ 


V^o-6 ’ \/6—^0 


, 0,1 


— § < wq < 0 


Ql3,14 (x c ,±-j&£ 


Ye V-ygoo-exg-, 




-v/ 6 —ujq 


,1 


■|<wo<0,^<^ 


Table 3. Existence conditions of the critical points at infinity of the system (2.41), located at the 

sphere X 2 + Y 2 + Z 2 = 1. We use the notation F 7 * 8 = ± 2 v / 3v / - 2 Acm)-3A L -+ -v^y-^ 0+3 ^ j g an 

v / A t /(2wo-9)(2wo+3)+5V4w 0 +6(3-2aj 0 ) 

arbitrary parameter. 


From the definition of z it follows that at the equilibrium point 


H - V(2'y-Xu)U 0 

V&y 


(2.65) 


Besides, under the conditions 7 > 0, Ay / 0, Xjj < 0 it follows that H necessarily tends to zero 
when the critical point is approached since <b tends to zero at the critical point. Furthermore, 
since the deceleration parameter satisfies q —>• — 1 at the critical point, it follows by continuity 
that H -C H 2 at late times, which means that H tends to zero too. 

Now, let us take advantage of the formula (2.51), which is valid up to fifth order, for 
obtaining some asymptotic expansions. Since = u 2 , then we get 

(A u ~ 27) 2 


7Af/Ay ln(a) + 2ci(Aj/ - 27)2 
Substituting (2.66) in (2.65) it follows 

^/TT~ ! 2 ~<- X U ( _ (Aq- 27) 2 _ 

N °V 7 V 7 -±Xy ln(a)+2ci(A(7-27) 2 


( 2 . 66 ) 


H = 


r-Af/ 

27 


\/6 


(2.67) 


Both expressions tend to zero as a ^ 00, provided 7 > 0, Ay / 0, Xu < 0. Integrating out 
(2.67) for a it follows 


a = exp 


(A u - 27 ) 2 2 a v-37 _ 


7(A[7-2 7 )^ 


2-^L-6c 2 


A[/A^(3 7 —A[/)(i 2 C/o(2 7 —Af/)—6 7 c|) 


2 2_ 

A t/-37 -|- 2ci 


7Ac/Av 


V 


~ e aitP1 (leading terms as t —>■ 00), 


( 2 . 68 ) 
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where c\ and C2 are integration constants, and 

2 x u -37 (Ay -2 1 ) x u-3~t +2 ( 12 - x u-'^ ( 3 Xu ~ 37 a v x u~'^ 


4t 1 

u'^~ Xu 


a\ = 


7-% 


and pi = ~ aF ?37 ’ w ^ ere a i > 0) and Pi > 0, provided 7 > 0, A// < 0. 

The expressions for the deceleration parameter is 

2t 1 I 7 (At/— 27) 2 (V6t*fU^J2 6c 2 ') i __2j_ 

2-^- -7) |1 


« = -! + ■ 


(Ac/ - 27)2 


(pi — lit - ? 1 

~ — 1-(leading terms as t —> 00). 

aiPi 


(2.69) 


Let us observe that a —>• 00 as t —>• 00. Additionally, it is recovered the expected effective 
equation of state parameter, u/ to t = — 1 as t —>• 00 for 7 > 0, Ay / 0, A// < 0. 

From the equations (2.39), (2.65), and (2.66), we obtain the following relation, which is 
valid at the critical point 


VW = 


TJ \ ( _ (Ay-27 ) 2 _ 

0 U V 7^u^v ln(“)+2ei(Ay— 27) 2 

27 


2-An 

7 


(2.70) 


Therefore, using the equations (2.66), (2.67), and (2.70), we can obtain asymptotic 
expressions for ( f>, H and </> in terms of t after the substitution of a ~ e Q1 * P1 , which are valid 
as t —> 00 for 7 > 0, Ay / 0, A u < 0. 

Now, in order for J4 to be a de Sitter solution it is required that pi = 1, which implies 
Xjj = 7. But this would lead to a contradiction, since 7 is assumed to be positive, and A u is 
negative in order for J4 to be an attractor. In conclusion, in the modified JBD theory, the de 
Sitter solution is not a natural attractor in the Jordan frame. Hence, we offer new arguments 
supporting the statements given in [26] against the validity of the results presented in [19-21], 
Furthermore, since J4 corresponds to T = 0, and since plays the role of an effective 
Planck mass, this would imply that there is no gravity at d> = 0. Thus, one must explicitly 
show that the solutions with a very small are suitable to describe the late-time universe. 
Indeed, we have 

R = 6H + 12 H 2 
= 6(1 — q)H 2 
„ / aa\ (a 

~ 12a 2 p 2 t 2 ( pi_1 ) + 6«i(pi — l)pit Pl-2 , (2-71) 

where in the last line we evaluated the approximate solution a — e aitP1 valid for 7 > 0, Ay 7^ 
0, A u < 0. On the other hand, the leading terms as a —>• 00 in (2.66) are 

(At/ — 2 7 ) 2 _ (At/-2 7 ) 2 __^ 

~ m 1 


$ -- 


'yXu^v ln ( a ) iXuXyaitP 1 
Henceforth, the term &R in the action (2.1) can be expressed as 


(2.72) 


$R ~ At pi ~ 2 + Bt ~ 2 , 


(2.73) 
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where A = - 12(At ^ A ^.“ lPl , B = - 6{Xu 1)pi , andpi = - 7^37 > 0, «i > 0, provided 

7 > 0, Xu < 0. Thus, &R tends to zero as the fixed point J 4 is approached, but it is always 
non-negative since A and B are both non-negative for 7 > 0, Xjj <0. Summarizing, we 
have explicitly shown that the solutions near J4 satisfy <f> ~ i, QR ~ At Pl ~ 2 + Bt~ 2 and 
H ~ aip\t pi ~ 1 , which tend asymptotically to zero but never reach this value. Besides, the 
effective gravitational coupling (the one measured in Cavendish-like experiments) [26], 


G e tt 


4 + 2a ; 0 ^_ 1 
3 T- 2cuo 


(2.74) 


satisfies 

G e s _ £ ^ Pi 

GeS t 


(2.75) 


As a consequence of the above, if we consider cosmological constraints on the variability of 
the gravitational constant [101], for instance the ones in [102], which uses WMAP-5yr data 
combined with SDSS power spectrum data: 


-1.75 x 10 -12 yr -1 < — < 1.05 x 10 _ 12 yr _1 , 

Cr 

or the ones derived in Ref. [103], where the dependence of the abundances of the D, 3 He, 
4 He, and ' Li upon the variation of G was analyzed: 

\G/G\ < 9 x 10 -13 yr -1 , 


one see that for a given p\ and a large t the above constraints are very easily fulfilled. Hence, 
there are solutions with very small <k which are suitable to describe the late-time universe. 

Now, let us use the above arguments to obtain values of the free parameters that lead 
to values of according to observations. Taking the Hubble time to be to = 13.817 x 10 9 
yr (as, for instance, in [20]), i.e., the present value of the Hubble constant Ho = 7.24 x 10 -11 
yr -1 , and given the value of A u < 0 , one gets 


18186.2 


ujq > —-2 -1.5 =7> -1.75 x 10 -12 yr -i < < 1.05 x 10 _lz yr“ 


-1 ^ G e g 


-,-12, -1 


and 


24911.1 , „ 

(jJ 0 > —-1.5 


A l 


G eS 


G e fi 


-13 ,„-l 


<9x10 16 yr 


Conversely, given ujq, 
190.716 


A u < — 


y/2uJo + 3 


-1.75 x 10 -12 yr -1 < ^ < 1.05 x 10” 12 yr ' 1 
C e ff 


and 


A[/ < 


223.209 
y/2bJo + 3 


G eff 


G e S 


< 9 x 10 -13 yr -1 . 
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2.4 The interplay with the induced gravity model 


Recall that our model (2.1) with U($) = Uo$,\u = 7 = —=,0 < 7 < i/f is equivalent 

V^o+f * 

to an extension of the so-called induced gravity model [76, 78] with action: 


Sjf = J V^g - ^-^-sTd^d^ - u 0 ® - \y v d^d v 4> - vw'J d 4 x. (2.76) 


Let us observe that J 4 exists for Xjj < 0 ,7 > 0 and the condition for recovering the 
induced gravity model is 7 = Xu, 0 < 7 < -y/|. Thus, strictly speaking, J 4 does not exist for 
the induced gravity model. 

In order to connect the main results of the section (2.2.1) with the results in the section 
( 2 . 1 . 1 ) we have to set V(cj>) = 0 in (2.76). Hence, the parameter space reduces to one 

parameter 0 < 7 < \l |, and the system (2.41) has the restriction 


x 2 + y 


1 1 , 


+ Z = 1. 


This restriction allows to eliminate one variable, say x. Additionally, the equation for e 
decouples. Thus, it is obtained the reduced dynamical system 


y' = \{y + 2 )y 2 - ^7 2 (y + 2 ) {{y + 2 ) 2 - 8 z 2 ) - 3yz 2 , 
z' = l (2 (y(y + 4) + 12)z - 3 7 2 (y + 2 fz) + 3 ( 7 2 - l) z 3 , 
defined on the invariant set 


(y,z),y 


1 1 , 


+ z 2 < 1, z > 0 > . 


(2.77a) 

(2.77b) 

(2.78) 


The only critical points that exist for (2.77) at the finite region are those enumerated below. 

1. J\ always exists. The eigenvalues of the linearization of the reduced system (2.77) are 

3 ^ ■ Thus, it is always the source. 

v6—37 v6—37 

2. J 2 always exists. The eigenvalues of the linearization of the reduced system (2.77) are 

■ Thus > it is alwa y s a saddle. 


3. The point J 7 becomes ( y,z ) = ( |. The eigenvalues of the linearization 


of the reduced system (2.77) are now -^|, ^ . Thus, it is the sink in the plane 


(y, z) for the range 0 < 7 < y 3 . 


2 4 


4 Although the point is unstable along the e-axis (then, a saddle) since t = ^ and y > 0 at the fixed point. 
However, the analysis on the e-axis can be dropped out since the corresponding equation decouples. 
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Concerning Jj, the results shown in the paragraph (2.2.1.5) cannot be applied since 
A u = 7 - Therefore, we introduce the new variable a = 2 -v / wo<J> to get 


2 a(t) 

VTTy 

y/3H(t )' 


(2.79a) 

(2.79b) 


Now, we want to demonstrated that the solutions given by (2.37) converges to Jj as 
t oo. 

First of all, since the relative errors (2.36) can be made small for large enough values of 
time, then we can approximate (2.37) by 


a(t) = e 
a(t ) = e 

Hit) = ' 


fz ( \ 7/67 __ 7/67 

v 2 7(C2 Ci > s inh 767+6 (A(t)) cosh 767-6 (A(t)), 

/3" / \ | - 7/67+2 7/67 — 2 

*^ C3 C2 C2 C3 sinh767+6 (A(r)) cosh 76^-6 (A(r)), 

■v/Cocsch(2A(t)) ( 27/67 — 3 ( 7 s — 2) cosh(2A(f))) 

3+2^6 - 7 2 


where A(t) = ^ / ^^ (24ci+t) . 

Secondly, substituting (2.80) in (2.79) we obtain 


67 (+6 — 7 cosh ( 2 A(f))) 

3 ( 7 2 — 2) cosh (2A(t)) — 27/67 

7 / 67/6 — 7 2 sinh (2A(t)) 

27/67 ~ 3 ( 7 2 — 2) cosh (2A(t)) 


2 7 2 


2 — 7 2 ’ 



2 — 7 2 


(2.80a) 

(2.80b) 

(2.80c) 


(2.81a) 

(2.81b) 


Thus, the point J 7 is approached by the solution (2.37) as the time goes forward. 

To end this section we substitute the values of y,z at J 7 in the equation (2.79). It 
follows by continuity that the solutions near J 7 satisfies the rates 




v{t)\j 7 = CT 0 e72V6-7 2 , 

(2.82a) 

( 2-7 2 )ty/U^ 


a{t)\j T =a 0 e TV 6 - 7 2 , 

(2.82b) 

HW u= ( 2 -rV?. 

7 / 27/6 - 7 2 

(2.82c) 

Now, comparing term by term the equations (2.82) and (2.80), and identifying <tq = 

e ^ 7(c2 - C3) ,ao 
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Cl = - 


it follows 


3 l/f 7 (c 3 -C 2 )+C 2 +C 3 


7 


In 4 6- ' 


V ( 6 —7 2 ){/q 


12V27 2 (7 0 


aftl - alt'll r -i2__V6 t+2 767-2 

--: —— 1 - = 1 — e v^Ve- 7 2 sinh v'i-y+e (A(f))cosh (A(t)) —> 0, 

a(t) 

cr(t) — a(t)\j 7 . . .. _*Sl. 


O-(f) 


= 1 — e V 2 V 6 - 1 2 sinhv'eT+e (A (t)) cosh ^-6 (A(t)) — >• 0 


if(t)-ff(i)|.7 7 = 


%/Uo (37 2 — 3 ( 7 s — 2) coth (A (t)) + 2 \/ 67 cscli (A (t)) — 6 ) 

3\/2 V 


0 . 


(2.83) 

(2.84) 

(2.85) 


Thus, we can use the asymptotics (2.82) as approximations for the solutions near J 7 . Sum¬ 
marizing, all the asymptotic results are consistent. 

In this example the “intermediate accelerated” solution does not exist, and the attractor 
solution has an asymptotic de Sitter-like evolution law for the scale factor. 


3 Field Equations in the Einstein’s frame 

Under the conformal transformation [104]: 

Qfiu = ^ 9 fiu where 4> = e 7X , (3.1) 

we can express the Jordan action (2.1) in the Einstein frame as: 

Sef = J v ^0 ~ - **(x) - ^e~ lx Q^d^(t)d v (t) - e^ 27 X U( 0 )^ d A x, (3.2) 

where il(x) = e _ 27 X I7(e 7X ) = Uoe~^ uX . In this frame x is interpreted as a conventional 
scalar field and 0 is now coupled with y. For a discussion of the equivalence between the two 
frames, see for example [105-116] and references therein. 

It is worth mentioning that in order for the conformal transformation to the Einstein’s 
frame be well defined it is required that the scalar field 4) satisfies 4> > 0. However, it can 
asymptotically evolves to its minimum value 4> = 0. This implies that x ~ 00 asymptot¬ 
ically, and also, as we will see in the next section, the dynamical system variable e tends 
asymptotically to zero at the fixed points. 

By considering a flat Friedmann-Lemaitre-Robertson-Walker (FLRW) metric 


ds 2 = — di 2 + a(t) 2 [dr 2 + r 2 ( d6 2 + sin 2 Odtp 2 )] , (3.3) 

the field equations become: 

X + Mx + |e~ 7X 0 2 ~ 2 7 e _27X U (0) + il'(x) = 0, (3.4a) 

<£ + 37*0 -70* + e- 7 *U'(0) = 0, (3.4b) 

37* 2 - Qx 2 + ^e _7X 0 2 + e- 27X U(0) + iX( X )) = 0, (3.4c) 

2 U + 37* 2 + Q * 2 + ^e _ 7 X 0 2 - il(x) - e- 27 X U(0)) = 0, (3.4d) 
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where now the dot means derivative with respect the conformal time t. 
The relations with the quantities in the Jordan frame is through: 


di 

dt 


V$, a = \/i>a, H = - 

a 


1 

\/¥ 


H + 


dt 




1 

\/¥ 


H + 


_1_d& 

2 $ dt 


We define the effective energy densities /^i, /^ 2 , the effective pressures pi and p 2 , and the 
coupling term Q as: [79]: 


di = \x 2 + U(x), (3.5a) 

fJ-2 = + e-^Vicj)), (3.5b) 

Pi = ^x 2 -U(x), (3.5c) 

p 2 = 1 e - e-^V^), (3.5d) 

Q = -'ye~ 2 ' rx x 4y(0) —e 7X 0 2 . (3.5e) 

Then, the system (3.4) is equivalent to 

fii + 3%(/ri + pi) = Q, (3.6a) 

/i 2 + %H{h 2 + P 2 ) = ~Q, (3.6b) 

T~i 2 = ^ (/L + ^ 2 ), (3.6c) 

Td = - 7 ; (di + Pi + d 2 + P 2 ), (3.6d) 


which is interpreted as the Einstein equations for two scalar fields coupled in a non-standard 
way. 


3.1 Dynamical system analysis 

In this section we present the above system by defining proper dynamical variables. Next, 
we examine the stability and discuss the properties of the solutions in this frame. 


3.1.1 Finite analysis 

Let us first note that the Friedmann equation (3.4c) can be expressed as: 

e-^0 2 , U(x) , X 2 , e-^V^) 


+ 


+ 


+ 


= 1 . 


m 2 m 2 m 2 m 2 

This expression suggests the introduction of the following variables 


1 e 2 x <j) 

e = e2 x , xi = 7=71 , yi = 


V6 n ’ 

The Friedmann equation (3.4c) leads to 

V(<t>) 


x _ V u (x) 

? z l 


Ven 


Vsu ' 


377 2 e 4 


1 2 , 2 , 2 1 

+ x l + y 1 +z( = 1. 


(3.7) 


(3.8) 


(3.9) 
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The evolution equations for the variables (3.8) are given by: 


1 = 7^1 (Sx\ + \/67?/i + 6 y\ -6) - Aye (xf + 2 /i + z\ - l) , 

y'\ = ^ ( 2 7 - 37^1 + \/6 (x? - 1) ?/i - 272/1 + %/62 /i + zf (At/ - 
4 = \zi ( 6 xj + Vi (%i - V / 6A//^ , 


e = V 2 ^ ie ’ 


which are defined on the phase space 

E := {(xi,yi,zi,e) € M 4 : x\ + 2 /? + z\ < 1 , 2:1 > 0 }, 

where the comma denotes derivative with respect f = in a. 

The relation with the variables (2.38) is given by 


e = e, xi = 

e = e, x = 


2x 

2 +? 

2 xi 

2 - \/ 672 /i 


2/1 = 


\y 


(2 + 22)7’ 

672/1 


zi = 


2z 

2 +V 


y = 


\/6 - 372/1 ’ 


2 = 


, with inverse 


22:1 


2 - -v/672/1 


(3.10a) 

(3.10b) 

(3.10c) 

(3.10d) 

(3.11) 


(3.12a) 

(3.12b) 


The conditions y 7 ^ 


—2 and 2/1 / 



guarantee that the direct transformation and its 


inverse given by Eqs. (3.12) are non singular. 

The relation between the coordinates defined in Einstein frame by (3.8) and the Poincare 
coordinates given by (2.61) is: 


X = — 


Y = 


2x\ 


( fa o\ /2(a:f+4v/672/l+2i-4) 

(v 6 7 2 /i - 2) yj \ yi{3iy ^ 2VE)+2 J 


+ 5 


672/1 


( rp. 0 \ / 2 (xf+4^72/1+^?-4) 

(V6 - 372/i) 


+ 5 


Z = — 


2z\ 


( fa. / 2 (xf+4+672/1+ 2 ?-4) 


+ 5 


The cosmological parameters are given by: 

A _ /+ _ 2 , 2 

“1 = TTTTo — 2/l + ^li 


m 2 

O — + 2 _ 1 2 2 

^ 2 = 

9 — 1 — -y2 = —1 + 3 (xf + 2/1)1 


ihtot = -1 - 


2 ?A 

37A 2 


= -1 + 2(x? + 2/i). 


(3.13a) 

(3.13b) 

(3.13c) 


(3.14a) 

(3.14b) 

(3.14c) 

(3.14d) 
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3.1.1.1 Fixed points and stability in the Einstein frame. The critical points, eigen¬ 
values and stability conditions for the system in the Einstein frame (3.10) are: 

1. E\ 2 : (x'i = 0,yi = ±l,£i = 0, e = 0). They always exist. The eigenvalues are: 


: \Z|t, ± \Z|t, 6 T 2\/67, 3 ^ y|A[/ 


£2 is a source for A u < \/ 6 , 0 < 7 < y §> or a saddle otherwise. £1 is always a saddle. 
These points are related with ,/i ,2 through (3.12). 


2. £3 : ^0, ^/| 7 ,0,0^. It exists for 0 < y 2 < §. The eigenvalues are: 

{7 2 ,2 7 2 -3,3 ( 7 2 -l),7(27-A r /)}. 

Whenever it exists it is always a saddle. This point is related to J 3 through (3.12). 

3. £4 : ^0,0, -A=^==, 0^. It exists for 7 > 0, Xjj < 0. The eigenvalues are 

{0, -3, ^ (-3 - v / 24 7 A[/ + 9) , ^ (-3 + y/24ry\u + 9 ) | . 

It is nonhyperbolic with a 3D stable manifold for A jj < 0,7 > 0. Thus, it has a large 
probability to attract the universe at late times. The full stability analysis requires the 
application of the center manifold theorem (the analysis is done in subsection 3.1.1.2). 
This point is analogous to J 4 studied before. 


4. £ 5 ^ : (±- 1 / 3 , 0 ,0,± —). They always exist. The eigenvalues are 


| — 1 , 2,1 (—Vl - 24 7 2 — l) , i (yi - 24 7 2 - l) | . 

Thus, they are always saddle. These points are completely analogous to J 5.6 examined 
before. 


5. £7 : ( 0 , ^=, y 1 — ^-,0 ). It exists for — y/E < A u < y/E. The eigenvalues 


are 


^(Ac /(7 + Ac/) - 6), - (Ac / 2 - 6) ,Xu(Xu ~ 27) j • 

Thus, it is a saddle. This point is related to J 7 through (3.12). 

6 . £ 8 , 9 : (±^^^,^ 7 ,^, 0 ). They exist for 0 < X v < y/E, 7 > ^ 
7 > 0, Ac/ > \/ 6 . The eigenvalues are 


or 


37 


A2 


A3 


7 + Ac/ ’ 2(7 + Ac/) 7 / 2 ’ 2(7 ± Ac/) 7 / 2 ’ 2(7 + Ac/) 7 / 2 J ’ 


- 31 




















where Ai .^,3 are the roots of: P( A) = I 27 A (( 7 s — 12 ) A u + 27 A fj + 67 + Xu 3 ) (7 + 
Xu) n/2 + AV 7 + ^t/ + 6 A 2 (57 - 2 Xu) (7 + A ^) 3 + 1447(27 - Xu) (Xu (7 + Xu) - 6)(7 + 
Xu) 9 - These points are related to Js,9 through (3.12). Observe that from the existence 
conditions it follows that the first eigenvalue is always positive, so, these critical points 
are either saddles or sources, but never an attractor. 

7. Eio,u ■ -, ~~, 0,0^ . They exist for — | < ujq < The eigenvalues are 

jl, ~ (l - a/ 25 - 24 7 2) , -i (l + a/25 - 24 7 2) , 2 - . 

They are always saddles. Since the transformation (2.38) is not well defined for yi = 
l~2 

A-4. there are no equivalent points in the Jordan frame at the finite region of the phase 
space (see the subsection 2.2.1.1). However, using the relation (3.13) we get that £ 10,11 
maps onto Q 3,4 which are located in the region at infinity (see Table 3). 

The existence and stability conditions for the critical points of (3.10) are displayed in 
Table 4. The relevant cosmological parameters and the description of the critical points in 
terms of their stability are shown in Table 5 5 * . 


5 As in the Jordan frame, the leading terms of the Hubble parameter in the neighborhood of the critical 

points are also shown. 
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Label 

xi,yi,zi,e 

Existence 

Stability 

El,2 

( 0 , ± 1 , 0 , 0 ) 

always 

Ei is a saddle 

£2 is a source for A u < y/6, wo > — |, 
saddle otherwise 

e 3 

( 0 ,y§ 7 , 0 , 0 ) 

0J0 > 

saddle 

E a 

f n n A/7 f± 

V ’ u ’ ±2 7 -a u ’ 7 

o 

VI 

b 

•< 

m|IN 

1 

A 

o 

3 

sink for 7 > 0 , Xy 7 ^ 0 , Xu < 0 . 

£ 5,6 

(±y/$,0,0,±%) 

always 

saddle 

e 7 


— \/6 < A u < \/6 

saddle 

00 

f + y/ X u(l+^u)-e ±6 ±7 n V ^ 3 

/l2wo+V24o;o+37+19 ^ \ ^ /H „ r 


^ 7 +Au ’7+A C/’VtTAE’^ 0> 2± 

\J 2oj 0 +3 S Af/ s VD or 

wo > — A{/ > \/6 

numerical inspection 

£io,n 

(=F'^ V f rI ,'f,0,oj 

3 ^ 1 

2 < Wq < 2 

saddle 


Table 4. The existence and stability conditions for the critical points of (3.10). 
























Label 



Q 

w tot 

Description 

m 

-El,2 

1 

0 

2 

1 

stiff matter 

Ho 

l+3Ho(t—to)' 

e 3 

27 2 

3 

1 272 

3 

2 7 2 - 1 

47 2 1 

3 

scaling solution. 

Ho 

l+ 27 2 -H 0 (t-to) ' 

E a 

2 7 

27 —Ay 

27-A(7 

-1 

-1 

Intermediate accelerated 

a(t) ~ e“ 2tP2 , «2 > 0 ,0 < P 2 < 1. 

~ a2P2^ >2 ~ 1 ■ 

Es,6 

0 

1 

1 

1 

3 

radiation-like. 

Ho 

l+ 2 H 0 (t-t 0 ) ' 

e 7 

1 

0 

i (A^ - 2 ) 

3 (^U ~ 3 ) 

quintessence-dominated. 

2Ho 

2 +A 2 ,Ho(t-to) ’ 

-E8,9 

7(7+Ai/)+6 

Ac/(7+Ac /)—6 

0 37 

1 27 

scaling solution. 

Ho 

(,1+^u)' 2 

(t+A uT 2 

“ 7+A(7 

7+A(7 

l+»(t-to)' 

Eio,ii 

2 

37 ^ 

1-2 

37 ^ 

1 

1 

3 

radiation-like. 

Ho 

l+ 2 H 0 (t—to)’ 


Table 5. Description of the cosmological parameters (3.14) of the critical points of (3.10). 





























3.1.1.2 Center manifold analysis for the intermediate accelerated solution E 4 

In order to investigate the stability of the center manifold for E 4 we introduce the new 
variables 


u = e, 


Vl = 


V2 = 


V 3 = 


XuXye 


+ xi, 


\/ 6 ( 2 7 - Xu) 
y/1 (V3 ~ VS^Xu + 3) 

V4 7 — 2Xjjy/&')Xu + 3 \/2 7 — A(7v / 8 7 A[/ + 3 

\/7 (V8jXu + 3 + \/3) 


+ 


VrAtm 


£1 (\/3 - yJ&'yXu + 3) 
2^8 jXu + 3 ’ 

"i (V8 7 A[/ + 3 + \/3) 


+ 


(3.15a) 

(3.15b) 

(3.15c) 

(3.15d) 


/V\ 

( 



«2 


V V 3 ) 

V 


\ 

/ U / 


(f\ 


Vl 

+ 

9i 


V 2 


92 

/ 

\v 3 ) 


\W 


(3.16) 


\/4 7 — ZXuVS'yXu + 3 v / 2 7 — Xu^/ 8 'yXu + 3 2y/8^Xu + 3 

to obtain 

/0 0 0 0 

0-3 0 0 

0 0 -i (3 + V24 7 Atf + 9) 0 

\0 0 0 -4(3-V24 7 A[/ + 9) ) 

where (/, <?i, <? 2 , 53 ) T is a vector of higher order terms. 

Since the center subspace of the origin is tangent to the e-axis, it follows that the center 
manifold of the origin is given locally by the graph 

{(a,ui,u 2 ,u 3 ) : vi = hi(u),v 2 = h 2 (u),v 3 = h 3 (u), 

Mo) = Mo) = Mo) = 0, 

h[( 0 ) = MO) = h’ 3 ( 0 ) = 0 , M < $}, (3.17) 

where 5 is a positive small enough constant. The functions hi,i = 1,2,3 satisfy a set of 
quasilinear ordinary differential equations which can be expressed symbolically as 

Kiuy-v't ] =0,i = 1,2,3, (3.18) 

L J Vi=hi(u) 

where one must substitute u', v \, v ' 2 , v' 3 through (3.15) and use the replacement v\ —> h\(u),v 2 - 
h 2 (u) and V 3 —>• h 3 {u). 

Setting 

hi(u) = a\\u 2 + ai 2 u 3 + O (u ) 4 , (3.19a) 

h 2 (u) = CL 21 U 2 + a 22 u 3 + 0 (a ) 4 , (3.19b) 

h 3 (w) = CJ 31 U 2 + a 32 u 3 + 0 (a ) 4 , (3.19c) 

in (3.18), equating to zero all the coefficients of equal powers of u, and solving for the o^-’s 

we get up to fourth order: 


an — 0 , ai 2 — — 


_ A[/ Xyj^yXu + 6 ) 

6 \/ 6 (A u - 2 7 ) 3 

^X 2 V X 2 V (V&yXu - VS'yXu + 3 + y/3) 


“ 21 \/2(2 7 - Xu) 5 / 2 (8\/3 7 A[/ + 3 (V8 7 A[/ + 3 + x/3)' 

(v / 3 7 A u + V87A u + 3 + \/3) 


031 = 


’ \/2(2 7 - Xu) 5 / 2 (8\/3 7 Atr - 3V8 7 Af/ + 3 + 3\/3) 



(3.20a) 

■,022 — 0 , 

(3.20b) 

P 

CO 

to 

II 

0 

(3.20c) 
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Henceforth, although the definition of the center manifold is quite different as for J4, the 
dynamics on the center manifold is given by the same equation (2.50), i.e., 


2 „,3 


u = 


'yX u XyW 


+ O {uf 


2(A u - 2 7 ) 2 

Neglecting the fifth-order terms and integrating we find that 


u(t) = ± 


\/-(Ac/ - 27 )" 


- 2 7 ) 2 + yAc/Ayf 


(3.21) 


(3.22) 


where ci is an integration constant that must be negative in order for u to be real-valued. 
Thus, as before, for 7 > 0, Xu ^ { 0 , 27 }, Avz / 0, it follows that the origin, and then £4, is 
stable provided Xu < 0 . 


3.1.1.3 Special case: A u = 0, !l(x) = Uq. 
new variables 

For A (7 = 0, i.e., it(y) = 

Uq we introduce the 

ui = e, 


(3.23a) 

U 2 = Zl~ 1, 


(3.23b) 

Vl = Xl, 


(3.23c) 

v 2 = yi + 2 ^ 

/I 7(21 " 

(3.23d) 

/ui \ /o 0 0 0 A 

u' 2 0 0 0 0 

v[ - 0 0-30 

\v' 2 J V 0 0 0 “3/ 

^ Ui \ 

U2 + higher order terms. 
Vl 

\v 2 ) 

(3.24) 


Since the center subspace of the origin is tangent to the plane ui-u 2 , it follows that the 
center manifold of the origin is given locally by the graph 


(3.25) 


u 2 ,v 1 ,v 2 ) : ui = hi(u 1 ,u 2 ),v 2 = h 2 (m,u 2 ), 
hr(0,0) = h 2 (0,0) = 0, Dh(0, 0) = 0, u\ + u\ < (5}, 


where Dh is the matrix of derivatives and 5 is a positive small enough constant. The func¬ 
tions h\, h 2 satisfy a set of quasilinear partial differential equations which can be expressed 
symbolically as 

' dhi(u ) , dhiiu ) , 


~Ui + 


-“2 - v i 


v i=hi(u) 


= 0 ,* = 1 , 2 , 


V dm 1 du 2 

where one must substitute u \, u ' 2 , v\, v ' 2 , through (3.24) and use the replacement v\ 
h\{ui,u 2 ),v 2 ->• h 2 (ui,u 2 ). 

Setting 


(3.26) 


hi = auu\ + ai 2 Uiu 2 + a 22 u% + 0(3), 
h 2 = bnu\ + bi 2 uiu 2 + b 22 u\ + 0(3), 


(3.27a) 

(3.27b) 
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and plugging back in (3.26), equating to zero all the coefficients of equal powers of u\ and 
U 2 , and solving for the a^-’s and ’s we get up to third order 


( 011 012 O22 
V &n bi 2 622 



Thus, the dynamics on the center manifold is dictated by 

u[ = -2^ 2 uiU2 + 0(3), 
«2 = 87 ^+ 0(3), 


(3.28) 


(3.29a) 

(3.29b) 


where 0(3) denotes error terms of third order in the vector norm. Neglecting the error terms 
and integrating out the system (3.29) we obtain 

ui = c 2 ^87 2 t + c 1 ,M2 = - 2 1 , -• (3.30) 

Finally, it follows that for A jj = 0, the origin, and then the point £4 behaves as a saddle 
since the orbits departs from the origin along the e-direction as the time goes forward. 

3.1.1.4 Features of the critical points of the system (3.10). Let us summarize the 
features of the critical points of the system (3.10) found in subsection 3.1.1.1: 


1 . Solutions E\ 2 are stiff-like solutions dominated by the kinetic term of the quintessence 
field x- This points are saddle so, they are not late-time solutions. 

2. £3 is a scaling solution where y 2 has the same order of magnitude than the effective 
potential U(<^>)/e 4 . 

3 . £4 is a solution where the potential energy il(y) of the quintessence scalar field and 
the effective potential U(c/>)/e 4 have the same order of magnitude. Besides wtot = — 1 , 
at the critical point. However, as we will prove later on subsection 3 . 2 , it does not 
represent generically a de Sitter solution. Indeed, the scale factor a satisfy 

a ~ e a2tP2 (leading terms as t —>■ 00), ( 3 . 31 ) 

where ct 2 > 0 and 0 < P2 < 1. 

4 . £5,6 is a saddle-like radiation dominated solution since w to t = |- 

5 . £7 represents the standard quintessence dominated solution. As a difference with the 
usual case [ 81 ], it can not attract the universe at late times since it is a saddle. 

6. £8,9 represents a scaling solutions where the energy density of both scalars scales with 
the same order of magnitud. Neither one field, nor the other, dominates the dynamics. 

7 . The critical points £10,11 mimics a radiation solution since tc?tot = v • Using the relation 
( 3 . 13 ) we get that £10,11 maps onto Q 3,4 (see Table 3 ). The existence conditions for 
these points lead to cao-values which are lower than the observed values for the BD 
parameter according to the recent analysis [ 22 - 25 ]. 
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3.1.2 Analysis at infinity 

Since the variable e is unbounded we add to aq, yi, z\ defined in (3.8) the new variable defined 
by E = A—, and the new time variable T given by dT = (1 — E)~ 1 dr. The resulting system 
admits the set of non-hyperbolic equilibrium points 6 at infinity x\ + y 2 + z 2 = 1 and E = 1. 
That is the boundary of the phase space H but with <f> —>• oo. 


3.2 Viability of the intermediate accelerated solution in the Einstein Frame 

In the Einstein frame, E 4 is probably the more interesting late-time solution. It represents 
a solution where the potential energy il(y) of the quintessence scalar field and the effective 
potential V(^)/e 4 have the same order of magnitude. Besides w tot = — 1 at the critical point. 
However, it does not represent generically a de Sitter solution in our scenario as we will prove 
in the following. 

From the definition of z\ it follows that at equilibrium point 


_ \J ( 2 q — Xu)Uq -AFy 


(3.32) 


•7 

Since e = e^ x tends to zero when £4 is approached and 7 > 0 , it follows that y a t 

late time. Now, under the conditions 7 > 0, Ay / 0, A u < 0, which imply the stability of £ 4 , 
it follows that T~L —>• 0 asymptotically. Furthermore, since the deceleration parameter satisfy 
q —>• —1 at the critical point, it follows by continuity that "H -C 7~L 2 at late time, which means 
that di tends to zero too. 

Now, let us take advantage of the formula (3.22), which is valid up to fifth order, for 
obtaining some asymptotic expansions. Since y = 7 Inu 2 , then we result in 


X(a) = - In 

7 


(Ac/-2 7 ) 2 

— 2 ci(A (7 - 27) 2 - 'yXuXy In a 


(3.33) 


Substituting (3.33) in (3.32) we obtain 


n 


t (2i - XuY 


2 -r (—2 ci (Xu - 27) 2 - 7 Xu Ay In a) 27 


(3.34) 


Integrating out for a it follows 

a = e / 32 +« 2 (r+t ) p 2 ^ e a 2 t p 2 ^ ea( jj n g terms as t —>• 00 ), ( 3 . 35 ) 


where 


a 2 = -- 


144 Ay-27 ( Ac7A 2^^_ 6 2 7 -a [/ ^(2 7 _A i7 ) 1 (u 0 (2-^-^ 2 t- a U 


2-, 

7F 5 ~ 


C ( 8 7 2 ^C/ I /— 

p 2 = — Xu x z 7 -r = v6e 2 (27-A u) J ^ an d ^ Notice that for Xu < 0,7 > 0 , 


£>( 2 -^ 


A(7-27‘ 


we have a 2 > 0,0 < p 2 < 1. This implies that a 00 as t —>■ 00 . 
This implies that the Hubble parameter is 

n = a 2P2 (T + i)^~ l 


(3.36) 


6 All the eigenvalues of the Jacobian matrix are zero. 
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and the deceleration parameter is 


(P2 - i)(r + t)- pz 

CH2P2 


(3.37) 


Additionally, it is recovered the expected effective equation of state parameter, wtot = — 1 as 
1 —>• 00 for 7 > 0, Ay / 0, Xjj < 0. 

From the equations (3.9), (3.32) and (3.33) we can obtain the following relation, which 
is valid at the critical point 


v(d>) 


tt \ ( _ (A;/ —27 ) 2 _ 

UpXu (2 7 -A;7)x _ ° U V 1^-u^y In(n)+2ci(Au—2 7 ) 2 

2y 27 


(3.38) 


Therefore, using the equations (3.33), (3.34), and (3.38), we can obtain asymptotic 
expressions for x,"H and <j) in terms of t after the substitution of a ~ e Q2tP2 , which are valid 
as t —>• 00 . For this point 


JK = 6 


1 + 


aa 


d 2 



2 


= na 2 2 p 2 2 (r + t) 2 ^-i) + 6a 2 {p 2 - i)p 2 (r + 


(3.39) 


and % = a 2 P 2 (r + t) P2 1 , where P 2 = ~aT? 2 7 ' Notice that for \jj < 0,7 > 0, we have 
07 > 0,0 < P 2 < 1. This implies that and T~L tends to zero as t —>• 00 in such a way that 
V\/U 2 -> 12 . 

Now, in order for £4 to be a de Sitter solution it is required that P 2 = 1, which implies 
A[/ = 0. But in this case the de Sitter solution will not be stable, rather, it will be a 
saddle, and the above approximation using the center manifold is not even valid. The other 
possibility for getting a de Sitter stage is to choose Xu < 0 and take the limit 7 —» + 00 . In 
this case the solution would be stable, but this would imply ujq —> — §. This value of cuo is 
several orders of magnitude lower than the bound wo > 4 x 10 4 imposed by the Solar System 
tests [22, 23], and the bounds estimated on the basis of cosmological arguments ujq > 120 
[24] and 10 < u 0 < 10 7 [25]. 


4 Concluding remarks 

In this paper we have investigated the case of a modified JBD theory which includes a 

2 _ 

power-law potential for the JBD scalar field explicitly given by t/(4>) = 7 , 7 = 

(uiq + 3/2)~ 4 ' 2 , where ujq is the BD parameter, and a quintessence field with an exponential 
potential V(4>) = as the matter content. This scenario was analyzed in the Jordan 

and Einstein frames and equivalences between them were discussed. We have presented the 
relation of our model with the induced gravity model with power-law potential and it was 
discussed the integrability of this kind of models when the additional quintessence field is 
massless, and has a small velocity. However, the main focus of this research was not to 
find analytical solutions as it was for the original induced gravity model, but the study of 
the asymptotic behavior of the solution space without using fine-tuning of the parameters 
and of the initial conditions. Nevertheless, the fact that for some fine-tuned values of the 
parameters we may get some integrable cosmological models, makes our choice of potentials 
very interesting. This issue deserves further investigation and it is left to future projects. 
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Secondly, using dynamical systems tools, we have provided the conditions on the pa¬ 
rameters of the theory that lead to an attractor with an effective equation of state parameter 
t^tot = — 1, describing the late-time evolution of the universe. The features of these solutions 
are essentially the same in both frames, i.e., analogous evolutions for the scale factor, and 
the conditions on the parameters are exactly the same. One condition is that the Brans- 
Dicke parameter satisfies wo > — §, which is compatible with the ranges described by the 
observations, and the complementary condition is that A u < 0 , i.e., a power-law potential 
U (4>) with a power greater than the second. These conditions are independent of the slope 
of the potential of the new scalar field, whenever its potential is not a constant. We have 
proved that de Sitter solution is not the natural attractor in the JBD model. In this model, 
the attractor condition is 7 = Xu, which is forbidden in our scenario. Instead, we have 
shown that the attractor in the Jordan frame corresponds to an intermediate solution of the 
form a(t) — e aitP1 as t —>• 00 where a\ > 0 and p\ = — y w ith 0 < p\ < 1 provided 
7 > 0, A jj <0. Furthermore, when we work in the Einstein frame we get that the attrac¬ 
tor is also an intermediate solution of the form a(t) ~ e" 2tP2 as t —>• 00 where 07 > 0 and 
p 2 = — y ~27 w dh 0 < P 2 < 1 for 7 > 0, Xu < 0. One possibility for getting a de Sitter stage 
is to choose Xu < 0 and take the limit 7 —>• + 00 . In this case the solution would be stable, 
but this would imply uo —>• — |, a value which is below of the lower limits on ojq imposed by 
the Solar System tests, and the bounds estimated on the basis of cosmological arguments. 

Furthermore, in the special case A u = 0, that is, for a quadratic potential in the Jordan 
frame, or for a constant potential in Einstein’s frame, the above intermediate power-law 
solutions are of a saddle type, which are not an attractor. These results were proved using 
the center manifold theorem, which is not based on linear approximation. 

Finally, for the specific elaboration of our extension of the induced gravity model in 
the Jordan frame, corresponding to the particular choice of a linear potential U(&) = FoT, 

Xu = 7 = (wo + 3/2) -1 / 2 ,0 < 7 < |, the dynamical system is then reduced to a two 

dimensional one, and the late-time attractor is linked with the exact solution found for 
the induced gravity model. This attractor solution satisfies the de Sitter-like asymptotic 

(2-7 2 )VCq 

expansion a ~ e V 2 V 6 - 7 2 as t —>• 00 , and the “intermediate accelerated” solution does not 
exist. 

Summarizing, apart from these special fine-tuned examples (the linear, and quadratic 
potential [/($)), it was shown that “intermediate accelerated” solutions are generic late-time 
attractors in our modified Jordan-Brans-Dicke theory. 
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